AO-AOB*  263 


NAVAL  SURFACE  WEAPONS  CENTER  DAHL OPEN  LAB  VA 
NSWC  SATRACK  MULTIPLE  SATELLITE  PROCESSOR  (SAMSAP). 


'  f 


, 


NSWC/DL-TR-3939 


CO 

<£> 

00 

00 


2ISKC  SATJBACK  MULTIPLE  SATELLITE 
PRO'CES SOR  ( SAMS AP )  , 


bu  . 

/A !  ^ .  JDWIN/iJOHNSON 
^Strategic  ~Svstem 


.strategic  Systems  Department 


/•J 


'  -// 


MARCH  1978 


Approved  for  public  release;  distribution  unlimited 


?>T 


'/■/''>  I  A-  •  i-t'j 


NAVAL  SURFACE  WEAPONS  CENTER 
Dahlgren,  Virginia  22448  / 


Silver  Spring,  Maryland  20910 


r  1 


1  r.  :  79 

OMB  NO.  22-R0334 


government-industry  data  exchange  program 

GENERAL  DOCUMENT  SUMMARY  SHEET 

Pt«os«  Ty pm  A(f  IflfonnofioB  —  S««  on  Rovors* 

1  OF  1 

1  ACCESS  NUMBER 

Elle^-o^ao 

2.  COMPONENT/PART  NAME  PER  GIDEP  SUBJECT  THESAURUS 

Computers,  Program,  Tape 

3  APPLICATION 

Ground 

«.  MFR  NOTIFICATION 

I  |  NOTIFIEO  [j{]  NOT  APPLICABLE 

5.  DOCUMENT  ISSUE  (Monr*/ V«0' 

March  1978 

6  ORIGINATOR'S  DOCUMENT  TITLE 

NSWC  SATRACK  Multiple  Satellite  Processor 
(SAMSAP) 

7.  DOCUMENT  TYPE 

fc~|  GEN  RPT  Q  NONSTD  PART  Q  SPEC 

B  ORIGINATOR  S  OOCUMENT  NUMBER 

TR-3939 

?.  ORIGINATOR'S  PART  NAME/lDENTl FlCATION 

N/A 

10  DOCUMENT  (SUPERSEDES*  (SUPPLEMENTS)  ACCESS  NO. 

None 

J.  ENVIRONMENTAL  EXPOSURE  CODES 

N/A 

12.  MANUFACTURER 

N/A 

3.  MANUFACTURER  PART  NUMBER 

. . . .  1 

1 15.  OUTLINE,  TABLE  OF  CONTENTS,  SUMMARY.  OR  EQUIVALENT  DESCRIPTION 

k 

Relates  to  performance  evaluation  of  fleet  ballistic  missile 
systems,  using  the  satellites  of  the  Global  Positioning  System 
for  error  analysis.  The  SAMSAP  computer  program  combines  the 
pass  normal  matrices  for  the  individual  satellites  into 
one  system  of  simultaneous  equations.  The  normal  equations 
may  be  expanded  to  include  the  parameters  of  the  station  and 
satellite  clock,  gravity,  and  polar  motion.  The  equations  are 
solved  for  the  parameter  corrections  and  the  systems  covariance 
matrix  at  epoch.  The  corrected  system  state  vector  and  the 
system  covariance  matrix  are  then  propagated  through  the 
trajectory  fit  span  and  output  to  a  tracked  ephemeris  data 
tape  that  is  used  by  Johns  Hopkins  University  in  their  Post- 
Flight  Processor. 

8  0  H  E 

i  101 

14  "BaTlis'tric  Missile;  Multiple  Satellite  Processor; 
System  Error  Analysis;  Ephemeris  Data 

Global  Positioning; 

(Doc  Des — M) 

17  GIDEP  REPRESENTATIVE  •• 

B.  A.  Howard,  Jr.  j 

P*«T!Cl%MtT  ACTIVITY  ^in  CQpfc.TT  —  -  _ 

U. S ►  Navy,  Naval  Sur¬ 
face  Weapons  Center,  Dahlgren  VA  (W4) 

00  2000 


REPRODUCTION  OR  DISPLAY  OF  THIS  '.-ATERIAL  FOR 
SALES  OR  PUBLICITY  PURPOSES  IS  PROHIBITED 


V 


: 

Lx 


INSTRUCTIONS  FOR  COMPLETING  THE  GIDEP  GENERAL  DOCUMENT  SUMMARY  SHEET 

NOTE.  Completion  of  a  Summary  Sheet  by  the  participant  is  not  mandatory  for  document  acceptance  into  GIDEP.  A  Summary  Sheet  will  be  prepared  by  the 
GIDEP  Operations  Center  for  document  submittals  received. 

PAGE  1  OF  Enter  the  total  number  of  summary  sheet  pages 
BOX 

1  Leave  blank — entry  will  be  completed  by  GIDEP  Operations  Center. 

2  Enter  standard  nomenclature  associated  with  GIDEP  Subject  Thesaurus  selected  from  Section  12,  Policies  and  Procedures  Manual. 

3  Indicate  application  which  the  device  was  used  (e.g.f  ground,  missile,  shipboard,  spacecraft;  refer  to  P  &  P  Manual,  Section  13). 

4  Device  manufacturer  must  be  notified  of  test  results.  Manufacturer  approval  of  report  is  not  required — include  pertinent  manufacturer 
correspondence  with  document  submittal  to  GIDEP;  check  NOTIFIED  entry.  If  document  is  for  a  nonstandard  part  or  of  a  general 
nature  and  a  manufacturer  is  not  identified,  check  NOT  APPLICABLE. 

5  Enter  month  and  year  of  document  issue. 

6  Enter  complete  document  title  exactly  as  it  oppears  on  originator  document. 

7  Identify  document  type  by  inserting  letter  X  by  appropriate  descriptor. 

8  Enter  document  number  exactly  as  it  appears  on  originator  document. 

9  Enter  part  name  and  identification  as  assigned  by  organization/agency  originating  the  report;  if  not  specified,  enter  N/A  (Not  Applicable). 

10  Delete  either  SUPERSEDES  or  SUPPLEMENTS.  If  document  supersedes/supplements  on  existing  GIDEP  document,  enter  GIDEP  microfilm 
access  number  of  appropriate  document.  If  document  neither  supersedes  nor  supplements  an  existing  GIDEP  document,  enter  the  word  NONE. 

11  Enter  the  single  symbol  coding  for  environmental  exposure  as  defined  in  Subject  Thesaurus,  Section  12,  P&P  Manual  (e.g.,  C  -  Salt  Spray, 

V  -  Vibration;  %  -  Shelf  Life);  if  not  specified,  enter  N/ A  (Not  Applicable). 

12  Enter  manufacturer  abbreviation  and  H-4  Code  number  listed  in  GIDEP  Manufacturer  List.  If  manufacturer  is  not  listed,  enter  the  phrase, 

SEE  BOX  15,  enter  manufacturer’s  full  name  and  division  (if  any)  in  Box  15.  If  more  than  one  manufacturer,  enter  phrase,  SEE  BOX  15,  enter 
additional  manufacturers  as  appropriate.  If  manufacturer  is  not  specified,  enter  N/A  (Not  Applicable). 

13  Enter  complete  part  number.  Use  open  0  for  alpha  letter  0,  and  use  0  for  numeric  zero.  If  more  thon  one  part  number,  enter  phrase,  SEE  BOX 

15;  enter  additional  port  number(s).  If  a  part  number  is  not  specified,  enter  four  dashes  ( - ). 

14  Enter  standard  part  number  such  as  the  IN  or  2N — diode  and  transistor  designators.  For  GIDEP  purposes,  any  military  ossigned  number  is 
considered  os  o  government  standard  part  number.  If  more  than  one  standard  number,  enter  phrase  SEE  BOX  15,  enter  additional  standard  number(s). 

15  If  subject  matter  in  document  can  be  categorized  into  more  than  one  major  subject  category,  enter  additional  sublet  categories  in  upper  right- 
hand  corner.  Briefly  summarize  test  results  or  material  detailed  in  text  of  document.  Include  any  pertinent  details  or  comments  required  for 
proper  interpretation  of  material  presented  (e.g.,  peculiar  environmental  capabilities,  unique  electrical  characteristics  that  may  be  "state-of-the- 
art,"  or  characteristics  thot  restrict  part  usage  to  particular  applications  or  any  other  details  that  may  aid  a  prospective  user  of  the  port). 

16  Enter  appropriate  words  or  phrases  that  enhance  information  retrieval  on  subject  matter(s)  contained  in  document.  As  a  secondary  data  retrieval 
technique  within  each  applicable  Major  Category  (entry  2),  the  document  is  referenced  in  the  computer  data  bank  and  Report  Index  according  to 
each  key  word.  Do  not  use  abbreviations  or  words  that  are  part  of  the  subject  category  listed  in  Box  2.  Key  word  phrases  are  limited  to  60  total 
characters  and  blonk  spaces.  Separate  key  words  and/or  phrases  with  commas. 

17  Enter  signature  or  name  of  GIDEP  Representative. 

18  Enter  name,  city,  and  state  of  participant  activity  or  corporation  and  division  submitting  the  document  and  GIDEP  two-character  code  (e.g.,  XI). 


DO  FORM  2000  (1  NOV  71)  (BACK) 


FOREWORD 
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INTRODUCTION 


This  report  presents  the  formulation  for  the  Naval  Surface  Weapons  Center  (NSWC) 
SATRACK  Multiple  Satellite  Processor  (SAMSAP).  SATRACK  is  a  Johns  Hopkins 
University/Applied  Physics  Laboratory  (JHU/APL)  program  to  evaluate  the  accuracy  of  U.S. 
Navy  fleet  ballistic  missile  systems.  The  program  uses  data  from  the  Global  Positioning 
System  (GPS)  satellites  in  conjunction  with  an  APL-developed  tracking  system.  The  satellite 
data  and  missile  telemetry  data  are  combined  to  obtain  an  accurate  estimation  of  the 
missile’s  trajectory. 

SAMSAP  is  an  orbit  determination  program  designed  to  meet  the  special  requirements 
for  generating  satellite  orbital  data  for  use  in  the  APL  Post-Flight  Processor  (PFP)  computer 
program.  SAMSAP  computes  a  multiple  satellite  solution  for  the  GPS  satellites  using  single 
satellite  pass  matricies  (range  and  range  difference)  and  perturbed  trajectories  from  the 
single  satellite  orbit  determination  program  CELEST.*  The  multiple  satellite  solution 
contains  solutions  for  station  and  satellite  clock  timing,  drift,  and  aging  errors.  It  may  also 
contain  solutions  for  up  to  20  gravity  parameters,  four  polar  motion  parameters,  and 
station  positions  (as  well  as  the  standard  orbit  constants).  Once  the  solution  is  obtained, 
the  ephemeris  for  each  satellite  and  the  system  covariance  matrix  are  propagated  as 
required  for  the  SATRACK  tracked  ephemeris  data  (TED)  tape.  These  results  are  written 
on  a  file  that  is  reformatted  and  output  on  magnetic  tape  for  transmittal  to  JHU/APL. 


The  CELEST  files  used  by  the  program  are 

1.  Pass  Matrix  Files  -  One  for  each  data  class  for  each  satellite.  This  file  contains 
the  pass  normal  equations  which  are  expanded  by  the  SAMSAP  program,  combined 
with  pass  matrices  of  other  satellites,  and  solved. 

2.  Perturbed  Trajectory  Files  —  One  for  each  satellite.  The  perturbed  trajectory  file 
contains  the  partial  derivatives  which  constitute  the  state  transition  matrix  for  the 
satellite.  The  partials  are  used  in  expanding  the  pass  matrices  for  gravity  parameters 
and  for  propagating  the  state  vector  and  covariance  matrix. 

3.  Polar  Motion  File  -  One  per  run.  This  file  contains  the  polar  motion  parameters 
needed  to  expand  the  pass  matrices  for  polar  motion. 

4.  ABCD  File  -  One  per  run.  This  file  contains  the  Greenwich  hour  angles  and  the 
inertial-to-earth  fixed  transformation  matrices  for  each  time  step  in  the  integration.  It 
is  used  by  the  SAMSAP  program  to  output  the  required  values  for  APL’s  use. 

‘James  W.  O’Toole.  CELEST  Computer  Program  for  Computing  Satellite  Orbits ,  Naval  Surface  Weapons  Center.  Dahlgren 
Laboratory  Technical  Report  NSWC/DL  TR-3565,  Dahlgren,  Virginia,  October  1976. 


Figure  1  is  a  flowchart  designed  to  illustrate  the  following  discussion  of  the  program’s 
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Figure  1.  Method  Flowchart 


The  program  takes  the  normal  equations  as  written  on  the  pass  matrix  files  and 
expands  them  to  include  station  and  satellite  timing.  If  required,  the  pass  matrices  can  be 
expanded  to  include  up  to  20  gravity  parameters  and  four  polar  motion  parameters.  The 
bias  section  of  the  matrices  are  then  eliminated. 

Two  different  tests  are  performed  to  determine  whether  or  not  the  matrices  should  be 
included  in  the  solution.  The  first  test  rotates  the  bias  section  of  the  pass  matrix  to  the 
radial,  tangential,  and  out-of-plane  directions  to  the  satellite’s  path.  The  rotated  equations 
are  solved  for  the  corrections  to  the  stations’  positions  (An,,  An2,  and  An3).  Each  of 
these  are  compared  to  an  input  tolerance  and  the  pass  is  rejected  if  any  one  of  the  tests 
fail.  The  second  test  consists  of  fitting,  in  a  least  squares  sense,  a  straight  line  through  the 
data  points  corresponding  to  An,  at  TCA  for  each  pass  in  each  rotated  direction.  The 
standard  deviation  of  the  fit  is  computed,  and  any  pass  which  lies  outside  of  some 
specified  multiple  of  this  standard  deviation  in  any  of  the  directions  tested  is  rejected. 

Ail  of  the  accepted  passes  are  then  summed  for  each  satellite.  The  dynamic  section  of 
the  accumulated  normal  equations  for  each  satellite  is  eliminated  and  the  resulting  equations 
are  summed. 

At  each  matrix  elimination,  the  reduced  variance  is  computed.  When  the  matrices  are 
summed,  the  variances  are  also  summed. 

After  all  matrices  have  been  accumulated,  the  resulting  normal  equations  are  solved  for 
the  parameter  corrections  (Ap’s).  The  Ap’s  and  covariance  corresponding  to  the  dynamics 
section  of  the  normal  equations  are  then  recovered.  These  Ap’s  are  used  to  adjust  the 
right-hand  side  of  the  bias  section  of  the  rotated  normal  equations  used  in  the  station 
navigation  section.  The  adjusted  equations  are  solved  for  the  correction  to  station  position 
(An,),  and  these  are  output  for  diagnostic  purposes.  Using  the  Ap’s  and  the  partials  from 
the  perturbed  trajectory  files,  the  satellites’  positions  are  propagated  over  the  fit  span  and 
output  to  the  TED  file  as  required.  The  covariance  is  also  propagated  over  the  fit  span 

and  written  to  the  TED  file  when  required.  For  each  time  to  which  the  covariance  is 

propagated,  the  satellite  position  section  of  the  covariance  is  rotated  to  a  reference  frame 
which  is  radial,  tangential,  and  out-of-plane  to  the  satellites’  trajectories  at  that  time.  The 
resulting  sigmas  in  the  rotated  frame,  as  well  as  the  sigmas  for  the  clock  terms  of  each 
satellite  and  station,  are  written  on  a  Diagnostic  Information  File. 

Since  the  TED  file  format  is  discussed  in  several  documents  as  well  as  this  report,  it 

will  not  be  detailed  here.  The  TED  file  written  by  this  program  must  be  reformatted  by  a 

second  program  into  a  form  usable  by  JHU/APL.  Appendix  A  contains  a  description  of  the 
final  tape  format. 
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SOLUTION 


Appendix  B  contains  a  list  of  parameters  which  should  be  included  as  input.  It  also 
contains  a  list  of  the  files  which  will  be  required  by  the  program.  Appendix  C  contains  a 
description  of  the  standard  CELEST  files  that  will  be  used  by  the  SAMSAP  program.  The 
logic  flowchart  (Figure  2)  has  been  numbered  to  correspond  to  the  numbered  steps  (1 
through  29)  in  this  report. 

On  the  pass  matrix  file,  the  parameters  will  be  in  the  following  order: 

Position 

Velocity 

Thrust  (if  present) 

Radiation  Pressure  (if  present) 

Station  Position 

Refraction 

Bias 

Drift 

Aging 

For  the  purpose  of  this  report,  position,  velocity,  thrust,  radiation  pressure,  and 
satellite  clock  parameters  (to  be  discussed  later)  will  be  referred  to  as  orbit  or  dynamic 
parameters.  Station  position  will  be  considered  both  a  bias  parameter  and  an  arc  parameter. 
Refraction  will  always  be  considered  a  bias  parameter.  Bias,  drift,  and  aging  (range  data 
only)  when  multiplied  by  c,  the  speed  of  light,  become  the  satellite  clock  parameters  and 
when  multiplied  by  minus  c,  become  the  station  clock  parameters.  Satellite  clock 
parameters  will  be  considered  dynamic  parameters,  and  station  clock  parameters  will  be 
considered  arc  parameters.  For  range  difference  data,  aging  is  not  present;  and  bias  and 
drift  refer  to  frequency  bias  and  frequency  drift,  respectively.  Frequency  bias  and  frequency 
drift  are  always  considered  bias  parameters. 


STEP  (1) 

Input  values  of  times  A6A,  A6B,  A5A,  and  A5B  (defined  in  Appendix  B)  are  packed 
as  follows: 


(YEAR  -  1900)  *  1000  +  DAY  +  SEC  *  10~5 


(1) 


Epoch  of  solution  is  A5A. 
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♦- 


Unpack  A5A,  A6A,  and  A6B  and  compute  seconds  from  A5A  to  A6A,  A6A  to  A6B, 
and  A5A  to  A6B. 

Let 


DAY5A  =  day  number  of  time  A5A 
YR5A  =  year  of  A5A 
SEC5A  =  seconds  of  time  ASA 

and  so  on  to  SEC6B. 

A5A  to  A6A:  If  YR6A  =  YR5A, 


SEC5T6  =  (seconds  from  A5A  to  A6A)  =  (DAY  6A  -  DAY  5A)  *  86400 

+  (SEC6A  -  SEC5A)  (2) 


If  YR6A  >  YR5A,  check  for  leap  year  [i.e.,  Mod  (1900  +  YR5A)  with  4];  if  it  is  0,  then 
YR5A  is  a  leap  year. 

For  leap  years. 


SEC5T6  =  (366  -  DAY5A  +  DAY6A)  *  86400  +  (SEC6A  -  SEC5A)  (3) 


For  non-leap  years. 


SEC5T6  =  (365  -  DAY5A  +  DAY6A)  *  86400  +  (SEC6A  -  SEC5A)  (4) 


Similarly  for  A6A  to  A6B.  ASA  to  A6B  is  then,  of  course,  the  sum  of  the  two 
preceding  calculations. 


STEP  (2) 

This  program  must  be  set  up  to  process  range  and/or  range  difference  data.  Although 
the  order  of  processing  is  up  to  the  programmer,  it  will  be  assumed  in  the  following 
discussion  that  the  order  will  be 

♦ 
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Range  data  for  first  satellite 

Range  difference  data  for  first  satellite 

Range  data  for  second  satellite  and  so  on 

Read  in  pass  from  i,h  satellite  from  pass  matiix  file.  Time-of-closest-approach  (TCA) 
on  the  file  is  given  in  seconds  from  the  beginning  of  the  day  in  record  type  2  of  the  file. 
Compute  the  seconds  from  ASA  to  TCA  (same  as  Step  (1)].  Check  to  see  if  this  is  within 
fhe  time  span  during  which  the  data -are  to  *be*  usei>  3f  not,  skip -to  next  pass.'  *•  •  ‘  *  * 


STEP  (3) 

Unpack  the  pass  matrix  and  form  single  satellite  normal  equations  (BAp  =  E). 
1 .  Range  data 


B  = 


Station  Position 
Station  Clock 
Dynamic 
♦Bias 


E  = 


Station  Position 
Station  Clock 
Dynamic 
♦Bias 


♦Bias  section  consists  of  refraction  and  station  position. 


Range  difference  data 


B  = 


Station  Position 
Dynamic 
♦Bias 


E  = 


Station  Position 
Dynamic 
♦Bias 


♦Bias  section  consists  of  station  position,  refraction,  and  frequency  terms. 
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At  this  point,  the  values  stored  in  the  station  clock  positions  should  be  the  negative 
of  the  terms  read  in  for  bias,  drift,  and  aging  (range  data).  The  values  stored  in  the 
satellite  clock  positions  should  equal  those  values  read  in  for  bias,  drift,  and  aging.  The 
suggested  order  for  arranging  dynamic  parameters  is 

Position 

Velocity 

Radiation  Pressure 
Thrust 
Clock  Bias 
Clock  Drift 
Clock  Aging 

STEP  (4) 

Multiply  all  rows  and  columns  of  B  and  rows  of  E  that  contain  clock  parameters  by 
c.  (Note:  Diagonal  elements  of  B  are  thus  multiplied  by  c2.) 


STEP  (5) 

Expand  the  normal  equations  for  gravity  parameters  specified  (from  0  to  20 
parameters).  Now  consider  the  normal  equations  partitioned  as 


where 

sub  a  =  arc  parameters 
sub  d  =  dynamic  parameters 
sub  b  =  bias  parameters 

Note  that  B  is  symmetric.  Thus, 

=  BT,  etc. 

a  a  da 
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Now,  letting  sub  g  indicate  gravity  parameters  and  d  indicate  the  first  six  dynamic 
parameters. 


aXTr*  T  3XTrA  ~T  _  aXTri  3XT 


r>  -'TCA  oaTCA  n  uaTCA  uaTCA 

B-  ’  :  ax—  aT 


ee 


ag  axn 


B  d  =  BT  =  T  ~T  B 

8d  d8  3g  axn 


dd 


R  -  RT  =  9^TCA  T  9XTCA  T  „ 

ga  "  ag  3g  3Xq  da 


(6) 


p  _  pT  _  aXTCA  T  3XTCA  T  n 

B-  '  Bl-  '  T  3X7“  ®db 


9gb  -  Dbg  3g 


_  axTCA  T  axTCA  ~T 
8  ag  axn  d 


where 


X0  =  epoch  elements 
superscript  T  =  transpose  of  matrix 
superscript  -T  =  transpose,  inverse  of  matrix 


dX-rcA  dXrcA 

3g  ’  ax0 


=  partials  obtained  from  the  perturbed  trajectory  file  for  the 
satellite;*  and,  X  represents  the  state  vector  of  the 
satellite 


* 

An  interpolation  tontine,  NTRP8,  mutt  be  used  to  obtain  all  values  from  the  perturbed  trajectory  file  (Appendix  D). 
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Thus,  the  normal  equations  can  now  be  written 


STEP  (6) 

Expand  normal  equations  for  polar  motion,  if  indicated. 
Compute 

r »  -«3 

A(TCA)  =  I  c«jj  1 

bco2  to, 


where 

to,  =  Aq  =  y 
<o2  =  Ap  =  x 
u>3  =  <o(At  +  tAt) 

At  =  UT1  -  UTC 

t  =  time  in  seconds  from  beginning  of  year 
to  =  rotational  rate  of  the  earth 

to  should  be  read  from  input.  Ap,  Aq,  and  At  are  available  for  the  start  of  each  day 
from  the  Polar  Position  File  (Appendix  E).  A  six-point  Lagrangian  interpolation  should  be 
used  to  obtain  Ap,  Aq,  and  At  at  TCA  of  the  pass.  At  should  currently  be  defaulted  to 
zero;  however,  provisions  should  be  made  to  read  it  from  the  Polar  Position  File  and 
interpolate  for  At(TCA)  using  the  six-point  Lagrangian  interpolator. 


Read  station  coordinates 


from  pass  matrix  file  and  compute 
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15 


STEP  (7) 


Now  consider  the  equations  partitioned  as 


"Bnn 

B„b“ 

APn” 

3“ 

_Bbn 

V 

-Pb_ 

3_ 

where  Bbb,  Apb,  and  Eb  are  defined  as  before  and  Bnn  contains  all  nonbias  terms. 


Compute  l/o2  for  bias  a’s  and  add  to  diagonal  terms  of  Bbb.  Compute  the 
bias-eliminated  normal  equations  and  reduced  variance: 


BEl  *  Bn„  -  BnbBb*Bbn 


EEl  =  En-BnbBb;Eb 


VEl  -  V-EjBb’Eb 


(12) 

(13) 

(14) 


Write  the  reduced  variance,  bias-eliminated  normal  equations,  and  TCA  of  the  pass  on  a 
bias  eliminated  file. 


STEP  (8) 

Perform  station  navigations.  Looking  at  the  bias  section  of  the  normal  equations 


Bbb^pb  Eb 


partition  the  matrices  as 


Bss 

B,c 

Es 

_b08 

Boo_ 

-Ap0_ 

3. 

13 
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where 

s  =  station  position  parameters 
o  =  all  other  parameters  in  the  bias  section 

Rotate  the  equations  to  the  radial,  along  track,  and  out-of-plane  directions  at  TCA.  The 
rotation  matrices  can  be  found  on  the  pass  matrix  files.  Let  r  be  the  3X3  rotation  matrix 
from  the  pass  matrix  files,  then  define 

R  ■  [o  °]  <161 

The  rotations  are  then 

Bbb  =  RTBbbR  O7) 

EbR  =  RTEb 

Now  compute 

An  =  (B^r'E*  (18) 

Output  An,  o’s  from  (Bj^)'1  (o’s  =  Vdiag  (B£b)_l ). 

Store  (B*br‘,  Eb  ,  R,  Bsp,  and  BS(J.  for  future  use. 

Bsp  =  station  position  -  polar  motion  submatrix 

Bsd.  =  station  position  -  dynamic  parameters  submatrix  without  the  station 
clock  terms 

STEP  (9) 

Perform  the  absolute  tolerance  test. 
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If 


An,  >  (absolute  tolerance), 

An2  >  (absolute  tolerance)2  (19) 

or 

An3  >  (absolute  tolerance)3 

tag  the  pass  bad  ((absolute  tolerance),  should  be  read  from  input).  Write  tag  (An,  2  3)  on 
bias  eliminated  file. 


STEP  (10) 

Return  to  Step  (2)  and  process  the  next  pass.  When  all  passes  from  the  i,h  satellite 
of  data  class  being  processed  have  been  completed,  go  to  Step  (11). 


STEP  (11) 

Perform  the  straight  line  test.  Read  from  the  bias-eliminated  file  TCA,  An,  2  3  for 
each  pass.  Form  the  following  matrices: 


M  = 

TCA>„..,  1  1 
TCAV.„  2  1 

N,  = 

^ni  Vass  1 
Ani>paSS  2 

JCA>p.„ .  !_ 

n2  - 

^n2^pass  1 

n3  = 

An3>pa„  r 

_A“2)pass 

A«Vp.» 

n 


Now  solve 


(MTMr‘(MTN,) 


=  (MTMr'(MTN2) 


(22) 


(MTMr1(MTN3) 


to  form  equations 


Let 


y  j  =  nij  t  +  b, 

y2  =  m2t  +  b2  (23) 

y3  =  m3t  +  b3 


*  =  TCA)pasS  j  j  *  1,2,  —  ,  n 

and  compute  y(,  y2,  and  y3  from  Equation  (23). 

Test 

IAn.  -  y(  .1  <  sig  mult  *  <r  j  =  1,  2,  •••  ,  N  (24) 

1  Pass  J  1  j=  i ,2,3 

where  sig  multj  is  read  from  input 
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and 


o.  =- 

N  An*  -  y,  2 

M2 

i 

_j=i  N 

i  =  1 ,2,3 


N  = 
sig  multj  = 


total  number  of  passes 
input  quantity 


(25) 


If  for  i  =  1,2,3  the  pass  fails  the  test,  tag  it  bad.  The  option  should  be  available  to  test 
on  i  =  1,2,  or  3  or  any  combination  thereof. 


STEP  (12) 

Accumulate  (paying  attention  to  station  and  timing  terms)  pass  matrices,  respecting 
tags,  for  the  satellite  being  processed.  Sum  elements  of  BE  t ,  Ee  i  ,  and  V1  •  as  indicated 
for  the  satellite  being  processed: 


BS  *  BE  1  (pass  |*BE|W  j-t  — ♦BEI)P, 

BS  -  EE|)p,„  N 

vS  -  vE')„„,  „ 


N 


(26) 


STEP  (13) 

Return  to  Step  (2)  and  begin  processing  next  data  class.  When  all  data  for  a  specific 
satellite  has  been  processed,  go  to  Step  (14). 


STEP  (14) 

Partition  the  summed  normal  equations  (Bs,  Es)  as 
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where 


d  -  dynamic  parameters 
m  =  all  nondynamic  parameters 

Add  l/o?  to  the  diagonals  of  the  dynamic  section  of  the  matrix.  Eliminate  the  dynamic 
section  of  the  equations. 


EEJ  -  E„  -En,dB;JEd 
VE:  =  Vs-  EXJEd 


Write  the  reduced  normal  equations  and  the  variance  on  a  dynamics  eliminated  file.  Save 
Bdd-  BddBdm’  and  BdiEd  for  future  use- 

Return  to  Step  (2)  and  begin  processing  data  for  next  satellite.  When  all  satellites  have 
been  processed,  go  to  Step  (15). 


STEP  (15) 

Accumulate  matrices  and  variances  as  before;  i.e., 

HA  -  0E2>.„  |+BEJ'„,  j  ♦-»' J>,lt  „ 

E*  -  EE=)„,  ,  +  •••  +EE!)„,  .  (29) 

VA  .  vEi)„,  ,  +  •••  „ 
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The  equations  are  now 


Polar  motion)A(j 

(if  present) 

^Ppolar  motion 

K 

polar  m otion 

1 

Polar  motion)At 
position)sta  , 

^Pstation  ptftftfcn 

mm 

F 

'  Station  posithffl 

-  •  ••  •  *«  ••  •  * 

i 

Position)sta  m 
Clock)sta  , 

^station  clock 

= 

V 

station  clock 

(30)  j 

Clock)4(a  m 
Gravity)  j 

^Pgravity 

E 

gravity 

•  (if  present) 

Gravity  )k 


STEP  (16) 

Add  in  the  appropriate  a’s  ( 1  /o? )  to  the  diagonal  terms  and  solve 


Apm  =  (BAr'EA  (31) 

Print  Apm  ,  (B A)~K 
STEP  (17) 

Recover  the  remainder  of  Ap’s  and  the  full  covariance: 


j  (APd>sat 

i 


<B;iEa 


sat  i 


(b:Jb 


d  d  dm 


>sat  iAPn 


for  satellite  i 


(32) 
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STEP  (18) 


Represent  the  full  covariance,  B" 1 ,  as 


Then, 


"K  Mt~ 
i  = 

_M  N  _ 


k  =  (BAr 


(33) 


M  = 


BddBdm 


'BddBdm  ^2 


BddBdm  \i 


:  n  =  number  of 
satellites 


(34) 


N  = 


Bdd\  0 

Bdi>2 


0  Bdi>n 


BddBdm  *1 
Bdd  Bdm  ^2 


BddBdm 


M1 


(35) 


Output  and  save  B  1 ,  and 


Ap  = 


Ap 

Apd) 


d'sat  1 


^Ph> 


d'sat  n 
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Compute  ami  print  total  reduced  variances  for  all  satellites  used  in  solution. 


(VE2)Pred .  =  EE2  .  -  Et2  Ap 

'  'sat  i  sat  i  sat  i  Fn 


(36A) 


(VA)Pred  =  yA  -  EA  Apm 


(36B) 


where 


Pred  =  predicted  values 


Let 


APn>sat 


'APm 

APd^sat  i 
APd\at  i 


STEP  (19) 

Perform  station  navigations.  Adjust  the  right-hand  side  of  the  rotated  equations 
[Step  (9) J  for  each  phase  of  each  satellite  by 


~b)adj 


= 


B 


sp  I  d 


APb 

APd' 


1 


sat  i 
pass  j 


where 


Bspid  "  f®sb  Bsd'f 

App  =  corrections  to  the  polar  motion  parameters 
Apd.  =  corrections  to  the  dynamic  parameters  (without  the  clock 
terms 
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Solve  for  Ani)jdj: 


Aiij) 


adj 


r'EbR) 


sat  i 
pass  j 


(38) 


Print  Anj)adj  . 

STEP  (20)) 

There  are  two  options  to  be  considered  at  this  point  that  depend  on  whether  or  not 
a  TED  file  is  to  be  written: 

1.  Propagating  the  satellites’  ephemerides  over  a  specified  interval  and  writing  the 
propagated  ephemerides  to  a  TED  over  a  portion  of  this  interval.  Propagating  the 
covariance  matrix  over  another  specified  interval  and  writing  the  covariance  to  the 
TED  file  at  specified  times  within  the  interval. 

2.  Propagating  the  ephemerides  and  covariance  over  the  specified  intervals  and  not 
writing  a  TED  file. 

In  order  to  propagate  the  ephemerides  and  covariance,  a  start  and  stop  time  and  step  size 
should  be  input  for  each.  In  addition,  a  start  and  stop  time  and  a  step  size  should  be 

input  for  the  interval  over  which  the  propagated  ephemerides  are  to  be  written  to  the  TED 

file.  It  can  be  assumed  that  this  interval  will  always  be  a  subinterval  of  the  time  span  over 
which  the  ephemerides  are  to  be  propagated.  The  step  size  can  be  assumed  to  be  an 
integer  multiple  of  the  step  size  used  in  propagating  the  ephemerides.  The  times  at  which 

the  covariance  matrix  is  to  be  written  to  the  TED  file  should  be  stored  in  an  array 

(TCOV,  Appendix  B),  and  can  be  assumed  to  always  be  within  the  interval  over  which  the 
covariance  is  being  propagated.  The  times  at  which  the  covariance  is  to  be  written  on  the 
TED  file  will  always  be  integer  multiples  of  the  step  size,  except  for  the  first  time  the 

covariance  is  to  be  written.  This  time  corresponds  to  missile  launch  time,  and  may  not  lie 

on  an  integer  multiple  of  the  step  size.  In  this  case,  a  separate  propagation  must  be 
performed.  At  each  time  line  of  the  covariance  propagation,  an  option  should  be  included 

to  write  a  diagnostic  reference  file  similar  to  Tape  73,  Appendix  C.  The  propagated 

ephemerides  and  the  covariance  matrices  written  to  the  TED  file  should  also  be  printed. 
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TED  FILE 


STEP  (21) 

The  first  record  written  on  the  TED  file  consists  of  two  words:  (1)  the  number  of 
satellites  being  processed,  and  (2)  the  number  of  time  lines  to  be  output  in  the  ephemeris 
section.  The  number  of  time  lines,  NTL,  is  given  by 


NTL  =  IFIX((A6B-  A6A)/output  interval]  +  1  (39) 


and  all  quantities  are  in  seconds. 

All  other  records  will  be  written  as  described  in  Appendix  F.  This  part  of  the  report 
will  cover  only  the  computations  required  to  obtain  the  desired  values. 

1  Header  record  -  no  computations  required 

1IA  Ephemeris  file  header  record  -  words  1  and  2  can  be  determined  from  input. 

Word  1  *  IFIX(A6A) 

Word  2  =  decimal  part  of  A6A  *  10s  (40) 


The  program  must  write  records  IIA  and  B  for  the  first  satellite,  then  the  second,  and 
so  on.  The  position  and  velocity  of  the  satellite  must  be  propagated  to  time  A6A  (see  the 
PROPAGATION  section).  The  right  ascension  of  Greenwich  and  the  transformation  matrix 
at  the  start  of  each  day  can  be  obtained  from  the  ABCD  file  (Appendix  G).  A  six-point 
Lagrangian  interpolation  should  be  used  to  determine  the  values  of  the  hour  angle  and 
transformation  matrix  at  the  time  needed.  Words  11  and  12  are  simply  read  from  input. 
For  the  first  satellite,  the  header  record  is  repeated  for  each  time  that  the  covariance 
matrix  is  to  be  output  in  part  IV  of  the  TED  format  [NCOV  times.  TCOV  (I)  contains 
the  times] .  The  only  words  that  must  be  changed  when  the  record  is  repeated  for 
satellite  1  are  words  10  and  13  through  21.  The  words  must  be  read  in  from  the  ABCD 
file  for  each  time  at  which  the  records  are  to  be  written. 


STEP  (22) 

Record  IIB  contains  satellite  position,  velocity,  and  clock  error  for  every  time  line 
between  A6A  and  A6B.  The  step  size  to  be  used  is  an  input  value.  Time  is  written  on  the 
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file  as  elapsed  time  in  seconds  from  A6A;  however,  all  calculations  must  be  done  in 
seconds  from  epoch,  A5A.  The  position  and  velocity  must  be  propagated  to  the  output 
time  and  the  clock  error  is  computed  by 


At  -  Apc)ock  bjas  +  Apclock  drift  ’  t  +  1/2  Apclock 


aging 


(41) 


where 


t  =  time  in  seconds  from  epoch 

Record  MB  is  repeated  for  each  time  line  until  A6B  is  reached.  Record  II A,  and  then  IIB 
are  written  for  the  next  satellite,  and  so  on.  If  records  II A  and  B  have  been  written  for 
all  satellites,  an  end-of-file  is  written. 


STEP  (23) 

Record  Type  III:  Station  Files  -  One  file  must  be  written  for  each  station,  and  the 
files  should  be  ordered  so  that  the  stations  occur  in  the  same  order  as  the  station  numbers 
are  input  (i.e.,  in  the  MSAT  array). 

All  words  contained  in  the  header  record  are  read  from  input  except  word  12  (the 
total  number  of  passes),  which  must  be  computed. 

To  write  the  station  files,  all  data  from  the  pass  matrix  files  must  be  sorted  by 
station,  accepted  or  rejected  based  on  TCA  of  each  pass,  and  counted.  This,  at  the  option 
of  the  programmer,  may  take  place  during  the  solution  part  of  the  program. 

From  the  pass  matrix  files,  which  are  satellite  ordered,  read  (for  each  pass)  year,  day, 
station  number,  and  number  of  points  in  the  pass  from  record  type  2;  start  (rise)  and 
stop  (set)  time  of  passes  from  words  five  and  six  of  record  type  4;  and  TCA  of  the  pass 
from  record  type  5.  As  in  the  solution  portion  of  this  program,  the  time  in  seconds  from 
A5A  to  TCA  of  the  pass  should  be  compared  to  the  time  spans  over  which  the  data  are 
to  be  processed.  Those  passes  which  lie  outside  the  interval  are  to  be  discarded.  For  those 
passes  which  are  accepted  TCA,  rise,  and  set  times  must  be  converted  to  seconds  from 
A6A  and  the  number  of  passes  accepted  for  each  station  must  be  counted.  When 
processing  range  or  range  difference  data,  all  passes  which  lie  in  the  acceptable  time  span 
are  to  be  counted.  When  processing  both  range  and  range  difference  data,  only  range 
difference  ^passes  are  to  be  tested  and  counted.  Output  would  then  be  a  header  record 
followed  by  the  appropriate  number  of  data  records  for  the  first  station  (an  end-of-file), 
and  then  the  header  record  for  the  second  station,  and  so  on. 
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STEP  (24) 


Record  Type  IV:  Solution  Parameters  and  Multiple  Satellite  Covariance  -  These  records 
are  written  for  each  time  specified  in  input  (TCOV  array). 

IV  A  Records  —  The  time  written  in  word  one  and  stored  in  the  TCOV  array  is  with 
respect  to  time  A6A.  To  propagate  the  state  vector,  it  is  necessary  to  reference  this  time 
to  epoch.  If  a  parameter  does  not  appear  in  the  satellites’  state  vector  (i.e. ,  radiation 
pressure  and/or  thrust),  default  values  (9.9999E9)  must  be  written  for  each  absent 
parameter. 

To  compute  the  appropriate  values  for  radiation  pressure  coefficients  and  thrust,  the 
initial  values  of  the  parameters  must  be  read  from  the  perturbed  trajectory  file  header 
record.  To  this  value  is  added  the  computed  correction  (Ap)  for  the  parameter;  i.e., 


K1(I)  =  K1(I)0  +  Apk  j(I) 
K2(I)  =  K2(I)0+Apk2(I) 
K3(I)  =  K3(I)0  +Apk3(l) 
T1  (I)  =  Tl(I)0+Ap„(I) 
T2(I)  =  T2(I)0  +  Apt2(I) 
T3(I)  =  T3(I)0  +  Ap(3(I) 


where 


K1(I),  K2(I),  K3(I)  =  x,  y,  and  z  components  of  the  radiation  pressure  coefficients, 
respectively,  for  the  i,h  satellite 

T1(I),  T2(I),  T3(I)  =  x,  y,  and  z  components  of  thrust,  respectively,  for  the  i,h 
satellite 

Subscript  0  =  initial  value,  as  read  from  the  perturbed  trajectory  file  of  the 
satellite  being  processed 

Apj(I)  =  computed  correction  for  the  jth  component  for  the  ith 
satellite. 

At  for  each  satellite  (station)  is  calculated  by  Equation  (41),  and  At,  At  for  each  satellite 
(station)  are  computed  from 


^Pclock  drift  +  ^Pclock 


aging 


t] 


sat 
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(sta  i) 


(42) 


£9 


At  -  APtlogfc  agjng  1  sa, 


(sta  i) 


The  a  priori  o’s  to  be  written  for  the  master  station  are  the  input  o’s  corresponding  to 
the  master  station  multiplied  by  1000. 


STEP  (25) 

Records  IVB  contain  the  station  positions  (x,  y,  and  z)  which  can  be  read  from  the 
pass  matrix  files.  The  refraction  correction  assumed  is  always  equal  to  zero.  The  a  priori 
o’s  in  words  3M  +  2  through  3M  +  4  are  the  input  a’s  multiplied  by  1000  which  correspond 
to  station  positions  x,  y,  and  z,  respectively,  in  the  bias  submatrix.  Word  3M  +  5  is  the 
input  a  for  refraction. 


STEP  (26) 

Records  IVC  contain  the  covariance  matrix  propagated  to  the  output  time.  The  matrix 
is  written  in  rows  for  the  upper  diagonal  half.  Each  write  record  contains  the  upper 
diagonal  half  of  a  row.  The  matrix  is  always  written  as  a  full  matrix;  i.e.,  it  has  locations 
for  all  parameters  whether  or  not  the  parameters  are  present  in  the  solution.  If  the 
parameters  are  not  in  the  solution,  default  values  are  written  in  all  rows  and  columns 
allocated  to  that  parameter. 


PROPAGATION 


STEP  (27) 

State  -  The  equation  for  propagating  the  state  is 


n  ax. 

tfi 


(43) 


where 


X,  =  corrected  values  at  time  t 

X,  =  uncorrected  values  at  time  t,  from  perturbed  trajectory  file 
n  =  total  number  of  parameters 
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Apj  is  obtained  from  the  solution  vector  of  the  satellite,  Apdj  in  Step  (17),  and  3X(/3pj( 

are  obtained  from  the  perturbed  trajectory  file.  If  thrust  is  present  in  the  solution,  the 
partials  of  the  orbit  parameters  with  respect  to  the  thrust  parameters  that  are  read  from 
the  perturbed  trajectory  files  must  be  adjusted. 


Let 


(t)  =  partials  of  the  six  orbit  parameters  with  respect  to  thrust,  as  read  from 
the  perturbed  trajectory  file  for  time  t 

\p(t)  =  partials  of  the  six  orbit  parameters  with  respect  to  position  and 
velocity  from  the  perturbed  trajectory  file  for  time  t 
ts  =  time  at  which  thrusting  starts,  from  perturbed  trajectory  file,  for 

satellite  i 

te  =  time  at  which  thrusting  ends,  from  perturbed  trajectory  file,  for 

satellite  i 

Then,  the  thrust  partials  required  for  propagation  of  the  state,  dX(/dT,o,  for  time  t  are 


3X, 

- —  =0  for  t  <  t 

3Tt0 


(43A) 


3Xt 

3T,„ 


^a(t)-  f(t)f~l(tsWa(ts)  for  ts  <  t  <  te 


(43B) 


T= —  =  Mt)W  l(teW/a(te)-  4>  for  t  >  te  (43C) 

d*t0 

It  should  be  noted  here  that  the  order  of  the  partial  derivatives  on  the  perturbed 
trajectory  files  does  not  correspond  to  the  ordering  suggested  earlier  in  this  report.  The 
suggested  ordering  corresponds  to  the  ordering  required  to  write  the  TED  files  and  was 
assumed  throughout  this  report.  However,  the  ordering  used  is  up  to  the  programmer.  If 
the  suggested  ordering  is  used,  care  must  be  taken  especially  when  propagating  the  state 
vector  and  the  covariance  matrix  to  insure  that  the  correct  partial  derivatives  are  used. 
State  in  this  case  refers  only  to  position  and  velocity. 


STEP  (28) 

To  propagate  the  covariance  matrix,  the  complete  state  transition  matrix,  i J/,  must  be 
formed.  \p  is  constructed  principally  from  the  partials  derivatives  found  on  the  perturbed 
trajectory  file  for  each  satellite.  The  equations  for  propagating  the  covariance  is,  in  general, 
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P(t)  =  «^(t)p(t0wT(t) 


(44) 


where 


P(tQ)  =  covariance  matrix  to  be  propagated 

In  this  instance,  however,  propagating  the  upper  part  of  the  covariance  matrix  containing 
polar  motion  and  station  position  is  not  desired.  Therefore,  P(tQ)  will  be  defined  as  the 
lower  square  partition  of  B~  *  when  it  is  partitioned  as 


where 


B'j1  =  station  position  portion  of  B~a 
The  station  clock  portion  of  Baa  is  contained  in  P(t0). 

\js  will  be  a  square  matrix.  It  can  be  constructed  as  follows: 


(45] 


sa  t  n 


(46A) 


^  = 


I  is  an  n  X  n  identity  matrix  where  n  is  the  number  of  gravity  parameters  present  in  the 
solution 


* 


sta  j 


1  t  1/2  t2 

0  1  t 

0  0  1 


(46B) 


t  =  time  from  epoch  to  which  covariance  is  being  propagated 
-  is  the  same  for  each  station. 


and 


*«i 


*u  = 


^sat 

i 

*  12 

0 

L_  *13 

axj 

3X, 

3X, 

3p, 

3p2 

•  •  •  PV 

3pL 

dX2 

3X2 

3X2 

3p, 

• 

3p2 

3pL 

• 

• 

• 

aXk 

• 

. 

3pj 

3Pl 

from 

the 

perturbed  i 

0 


(47) 


(48) 


\pn  may  or  may  not  be  square,  depending  on  L  (the  number  of  parameters  in  the 
solution).  If  thrust  is  a  parameter,  the  partials  with  respect  to  thrust  must  be  adjusted  in 
the  manner  described  in  Step  (27),  Equations  (43A),  (43B),  and  (43C). 


^/j2  =  1,  the  identity  matrix 


(49) 
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The  dimensions  of  \p{2  are  dependent  on  whether  or  not  thrust  and/or  radiation  pressure 
parameters  are  present  in  the  solution.  If  thrust  or  radiation  pressure  alone  is  present,  i|/j2 
is  a  3  X  3  matrix;  if  both  are  present,  i//j2  is  a  6  X  6  matrix 

^  ,  .  (50) 

*i3  *sta  J 

^j3,  like  i pst3  j,  is  the  same  for  all  satellites. 


STEP  (29) 

Covariance  Statistics  —  To  compute  the  covariance  statistics,  the  state  of  each  satellite 
must  be  propagated  to  the  time  at  which  the  statistics  are  to  be  computed  [Step  (27)]. 

A  A  A 

The  rotation  matrices  Hj  =  (R,  I\  A);  can  be  formed  from  the  state  vector 
where 


V.  =  (Rt  X  R.)  X  Rj 


A.  =  Ri  X  R, 

RM  =  Xsa «  i 

RXi 

=  *sa.  i 

Ry  1  ”  ^sat  i 

Ry. 

=  ysa.  i 

~  Zsat  i 

K 

=  Zsat  i 

A  R; 

=  (R2  +  R2  +  R2  )1/2 
xi  yj  z i 

Ri  =  W~\ 

lR,l 

r 5- 

■  ir.i 

lf,l 

=  (r2  +r2  +  r2  ),/2 

^ i  y  i  c\ 

A  Aj 

Aj  =  K\ 

|Ajl 

=  (A2  +  A2  +  A2  )1/2 
x,  Vj  Zj 

i  =  1,2,  •••  ,n 
n  =  number  of  satellites 
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Using  the  rotation  matrices,  rotate  the  position-position  part  of  the  propagated  covariance 
matrix  for  each  satellite. 


pK(t)  =  HjPpfOHj 


where 


Pp.  =  position-position  part  of  P(t)  for  the  ilh  satellite 

The  rotated  values  should  be  printed.  The  square  root  of  the  diagonal  terms  of  Pp  and 
the  square  roots  of  the  diagonal  terms  corresponding  to  the  clock  bias,  drift,  and  aging  of 
the  satellites  and  the  stations  should  be  written  on  the  Diagnostic  Information  File 
(Appendix  H). 
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APPENDIX  A 

SATRACK  FINAL  TAPE  FORMAT 
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MEMORANDUM 

From;  NSWC/DL  (E.  Johnson) 

To:  APL  (R,  J,  McConahy) 

Date:  25  January  1977 

Subj:  Comments  on  "SATRACK  File  Definition  Sheet" 

Ref;  (a)  "SATRACK  File  Definition  Sheet,"  Version  1,  dtd  11/18/76 

(b)  "SATRACK  Ephemeris  Data  Tape  Binary  Format,"  1/20/77  Revision 

1.  The  "SATRACK  File  Definition  Sheet"  (Ref.  (a))  as  furnished  to  NSWC/DL 
by  APL  exhibits  several  shortcomings  which  quite  probably  result  from  in¬ 
complete  knowledge  of  the  logic  in  the  TED  tape  generating  programs  at 
NSWC/DL.  It  was  thus  felt  that  it  would  be  better  to  thoroughly  document 
this  logic  for  APL's  use,  rather  than  commenting  on  reference  (a)  line-by¬ 
line. 

Those  items  in  reference  (a)  which  appear  to  be  correct  and  should 
not  be  affected  by  any  required  changes  in  APL's  program  are: 

(a)  the  listing  of  the  possible  file  sequences,  and 

(b)  the  formats  of  those  records  which  are  of  fixed  length. 

In  the  following  discussion  of  the  tape  format,  tape  records  will  be 
referred  to  as  denoted  in  reference  (b).  A  diagram  of  the  tape  format 
has  been  included  to  further  clarify  the  tape  structure. 

Packing  words  used  in  records  I-IIIB  are  either  6's  or  7's.  The  6's 
are  used  only  for  header  records,  and  the  7's  are  used  for  data  records. 

In  part  IV  the  packing  word  takes  the  form  747X7YY7ZZZZ  where  X  =  1 ,  2  or  3 
Indicating  record  IVA,  IVB  or  IVC.  respectively;  YY  =  01,  02...,  nCOV  in¬ 
dicates  the  number  of  times  record  type  IV  has  been  repeated;  and  ZZZZ  =  3010 
for  parts  IVA  and  IVB.  In  part  IVC  the  covariance  matrix  is  written  row  wise 
for  the  upper  diagonal  half,  and  ZZZZ'  indicates  that  the  values  that  follow 
comprise  that  row  of  the  upper  diagonal  half  of  the  covariance  matrix  which 
contains  ZZZZ  number  of  elements.  Since  the  covariance  matrix  is  square, 

ZZZZ  is  also  equal  to  the  total  number  of  rows  in  the  covariance  matrix  +1 
minus  the  row  number  being  written,  i.e.,  15N+3M+l-i. 
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To  facilitate  further  discussion  of  the  tape  format,  it  is  convenient 
to  define  a  "packing  record."  A  packing  record  consists  of  a  packing 
word  and  nine  data  or  padding  words.  If  the  data  being  printed  does  not 
completely  fill  the  nine  allocated  words,  the  remaining  words  are  filled 
with  padding  words  ( .99991 00000000E+10).  All  data  on  the  tape  is  con¬ 
tained  in  packing  records.  It  should  be  noted  that  the  nine  data  words 
are  counted  as  per  reference  (b);  thus,  the  fact  that  all  alphanumeric 
words  in  reference  (b)  are  two  element  arrays  and  appear  on  the  TED  tape 
as  a  single  word  must  be  accounted  for. 

The  tape  is  written  in  blocks  of  200  words.  The  blocks,  however, 
cannot  be  interpreted  strictly  as  logical  records  or  as  representing  a 
particular  portion  of  the  TED  tape.  If  the  program  writing  the  TED  tape 
encounters  an  end-of-life  (as  specified  in  reference  (b))  and  the  data 
being  written  has  not  filled  a  200  word  block,  the  remainder  of  the  block 
is  padded  out. 

The  header  record  (Part  I)  of  the  TED  tape  is  a  variable  length 
record  depending  only  on  the  number  of  satellites.  Packing  records  are 
written  until  all  of  the  data  in  Part  I  is  written.  If  the  last  data 
word  in  Part  I  does  not  correspond  to  the  10th  word  of  a  packing  record, 
the  packing  record  is  padded  out. 

The  next  word  written  is  the  packing  word  for  the  1st  packing  record 
of  Part  IIA.  This  is  a  constant  length  record  of  21  words;  thus,  three 
packing  records  are  written,  the  last  of  which  is  padded  out.  For 
satellite  #1,  Part  IIA  is  repeated  NC0V  times.  The  packing  record  follow¬ 
ing  the. (NC0V+1 )  IIA  record  begins  with  a  7  and  contains  the  eight  data 
words  of  I  IB  and  one  padding  word.  These  records  are  repeated  for  each' 
time  line  in  the  ephemeris  output  span.  Record  IIA  is  then  written  (only 
once)  for  the  2nd  satellite;  records  I  IB  are  written  for  the  2nd  satellite 
and  so  onto  satellite  N.  If  the  last  I IB  packing  record  written  for 
satellite  N  does  not  correspond  to  the  end  of  a  200  word  block,  the  re¬ 
mainder  cf  the  block  is  padded  out. 

Record  I  IIA  starts  at  the  beginning  of  the  next  200  word  block.  It  is 
a  fixed  length  record  of  12  words;  therefore,  it  always  appears  in  two 
packing  records.  The  last  six  words  of  the  second  packing  record  are 
padded  out. 

The  next  packing  record  contains  the  6  data  words  of  record  1 1 1 B  and 
3  padding  words.  These  records  are  repeated  until  all  passes  have  been 
listed.  If  the  last  1 1  IB  record  written  does  not  correspond  to  the  end  of 
a  200  word  block,  the  block  is  padded  out. 

Records  I IIA  ard  1 1 IB  as  described  above  are  then  written  for  station  #2 
and  so  on  to  station  M. 
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Record  IVA  starts  at  the  beginning  of  the  next  200  word  block,  the 
first  packing  word  being  written  as  747170170000.  Since  this  is  a  var¬ 
iable  size  record  of  considerable  length  depending  on  the  number  of  sat¬ 
ellites  and  the  number  of  stations,  several  packing  records  are  required 
to  completely  write  the  record  on  the  tape.  The  packing  word  in  all  of 
these  packing  records  is  747170170000  indicating  record  IVA,  1st  time. 

When  the  last  packing  record  of  IVA  is  written  and  padded  out  as  required, 
the  packing  word  is  changed  to  747270170000  indicating  record  IVB,  1st 
time.  This  is  also  a  variable  length  record  requiring  a  variable  number 
of  packing  records;  the  last  of  which  may  be  padded  out. 

After  the  last  packing  record  for  part  IVB,  the  packing  word  is 
changed  to  74737017XXXX,  indicating  part  IVC,  1st  time,  and  XXXX  as  in 
the  preceding  discussion  of  packing  words.  Row  1  of  the  covariance  matrix 
is  then  written  using  the  number  of  packing  records  required  by  the  size 
of  the  matrix.  The  packing  word  is  then  changed  to  74737017XXXY,  where 
XXXY  =  XXXX  -1,  which  indicates  that  the  row  of  the  upper  diagonal  half  of 
the  covariance  matrix  now  being  written  contains  1  less  element  than  the 
preceding  row.  The  number  of  packing  records  used  for  each  row  can  be 
computed  from  the  value  XXXX  in  the  packing  word  for  that  row.  The  process 
is  repeated  until  the  packing  word  is  747370170001  indicating  that  the 
last  row  of  the  covariance  matrix  has  been  written.  Part  IVA  is  then 
written  for  the  2nd  covariance  output  time,  and  the  packing  word  for  this 
cycle  is  747170270000. 

Parts  IVB  and  IVC  are  then  written,  as  described  above,  for  the  second 
covariance  output  time.  This  process  is  continued  until  the  packing  word 
becomes  74737XX70001 ,  indicating  the  last  row  of  the  covariance  has  been 
written  and  XX  =  nCOV  indicates  that  the  covariance  has  been  written  for 
all  the  time  lines  requested.  If  the  last  packing  record  written  does 
not  correspond  to  the  end  of  a  200  word  block,  the  remainder  of  the  block 
is  padded  out. 

There  are  no  locations  on  the  tape  in  which  default  words 
(.9999900000000E10)  are  used  exclusively.  Default  words  fill  spaces  that 
are  assigned  to  quantities  that  have  not  been  used  or  solved  for  in  the 
particular  run,  but  may  appear  in  other  runs. 
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TED  TAPE  FILE  DIAGRAM 
(Based  on  Reference  (b)) 


1.  Tape  Header  I 

Satellite  1  H 


Satellite  2  11 


Tape  Header 


A 

A 


|  NCOV  +  Times 


A 


B 

B 

One  B  for  Each  Time  Line 


B 

A 


B 

B 


One  B  for  Each  Time  Line 


Satellite  N 


Station  I 

i 

i 


ii 


in 


B 

B 


B 

E-O-F 

A 

B 
B 


One  B  for  Each  Time  Line 


One  B  for  each  Pass 
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TED  TAPE  FILE  DIAGRAM  (Continued) 

(Based  on  Reference  (b)) 

E-O-F 


III 


IV 


IV 


A 


B 

B 


One  B  for  Each  Pass 


B 


i 


E-O-F 


A 

B 

C 


A 

B 

C 


NCOV  Times 


E-O-F 

E-O-F 
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TRACKED  EPHEMERIS  DATA  TAPE 
BINARY  FORMAT 


Revised  1/19/77 


Tape  Header  (1  Per  Tape) 


Word  Number  Format 


1,2  A 

3,4  A 

5,6  A 

7  I 


8 

8+  1 
8  +  2 


8  +  N  1 

8  +  N  +  1  I 

8  +  N  +  2,  8  +  N  +  3  A 

8  +  N  +  4  1 

8  +  N  +  5  E 

8  +  N+6  E 

8+N+7  E 

8+N+8  E 


Definition 


APL  -  Tape  User 
NSWC/Dahlgren  -  Tape  Generator 
MODX  -  OD/CT  Version  X 
(Yr  -  1900)  *  1000  +  Day  Number 
Run  Specified  (Al  -  RSF) 

N  -  Number  Satellites 
N1  -  Satellite  Number 
N7  -  Satellite  Number 


NN  -  Satellite  Number 
M  -  Number  Stations 
RW  or  FW  -  Real  or  Fitted 
Trajectories 

NCOV  -  Number  of  Covariance 
Matrices  Written  Out 
(Yr  -  1900)  *  1000  +  Day  Number 
+  Sec/105  (A5a-RSF) 

(Yr  -  1900)  *  1000  +  Day  Number 
+  Sec/105  (A5b-RSF) 

(Yr  -  1900)  *  1000+  Day  Number 
+  Sec/  10s  (A6a-RSF) 

(Yr  -  1900)  *  1000  +  Day  Number 
+  Sec/105  (A6b-RSF) 


II. 


Ephemeris  Files  (1  Per  Satellite) 


A.  Header  Record*  1  * 


Word  Number 

Format 

Definition 

1 

1 

(Yr-  1900)  *  1000+  Day  of 

Ephemeris  Record 

2 

E 

t  -  Seconds  UTC  After  Midnite 
(A6a-RSF ) 

3 

I 

Satellite 

Number 

4 

E 

X  Satellite  Number  (A6a-RSF) 

5 

E 

Y  Satellite  Number 

6 

E 

Z  Satellite  Number 

7 

E 

X  Satellite  Number 

8 

E 

Y  Satellite  Number 

9 

E 

Z  Satellite  Number 

10 

E 

True  Right  Ascension  of  Greenwich 
at  Epoch  A6A-PSF  (Radians) 

11 

E 

X  —  Complete  Polar  Motion  at  Epoch 

12 

E 

Y  -  Complete  Polar  Motion  at  Epoch 

13 

E 

(1,  1)  \ 

i 

14 

E 

(1,  2) 

15 

E 

(1,  3)  , 

1  Elements  ot  Transformation 

16 

E 

(2,  1) 

Matrix  From  Mean  1950.0 

17 

E 

(2,  2) 

>  Inertial  to  Earth -Fixed 

18 

E 

(2,  3) 

Geocentric  Equatorial  at 

19 

E 

(3,  1) 

20 

E 

(3,  2) 

|  Solution  Epoch  (A6a-RSF) 

21 

E 

(3.  3)  ) 

^Note;  Header  record  written  for  time  A6a  for  all  satellites;  for  satellite  number  1,  it  is  repeated  NCOV  times  with  the  tight 
ascension  of  Greenwich  and  the  transformation  matrix  updated  to  correspond  to  the  times  that  the  covariance  matrix  will  be  output 
in  part  IV. 
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B.  Ephemeris  Data  Records  (1  Per  Time  Line) 


Word  Number 

Format 

Definition 

1 

E 

Time  (sec)  Relative  to 
Epoch  (A6A) 

2 

E 

X  —  Coordinate 
Mean  1950.0 

(km)  in 

3 

E 

Y  -  Coordinate 
Mean  1950.0 

(km)  in 

4 

E 

Z  —  Coordinate 
Mean  1950.0 

(km)  in 

5 

E 

X  -  Coordinate 

Mean  1950.0 

(km)  in 

6 

E 

Y  -  Coordinate 
Mean  1950.0 

(km)  in 

7 

E 

Z  -  Coordinate 
Mean  1950.0 

(km)  in 

8 

E 

At  -  Satellite  Clock  Error 

The  last  epheremis  file  terminates  with  E.O.F. 
UI.  Station  Files  (1  Per  Station) 

A.  Header  Record 

Word  Number  Format 


1 

2 


I 

E 


3 


E 


4 

5 


E 

1 


6 

7 


E 

E 


8 


E 


Definition 


Station  Number 
<t>  —  Latitude  (deg)  in 
WGS-72  System 
X  —  Longitude  (deg)  in 
WGS-72  System 

h  -  Height  (m )  in  WGS-72  System 
Data  Class 

1  =  Range 

2  =  Range  Difference 

3  =  Both 

Mean  on  Range  Data  (=  0.0) 
Sigma  of  Noise  on  Range  Data 
in  Data  Generation  R 
Range  Data  Rate  -  Data 
Generation  Run  (sec) 
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Station  Files  (1  Per  Station)  (Continued) 


Word  Number 

Format 

Definition 

9 

E 

Mean  on  Range  Difference  Data 
(=  0.0) 

10 

E 

Sigma  of  Noise  on  Range  Difference 
Data  in  Data  Generation  R 

11 

E 

Range  Difference  Data  Rate  - 
Data  Generation  Run 

12 

1 

K  —  Total  Passes  From  Station 

B.  Data  Records  (1 

Per  Pass) 

1 

1 

Pass  Number 

2 

E 

Start  Time  (sec  from  A6a) 

3 

E 

End  Time  (sec  from  A6a) 

4 

I 

Tracked  Satellite’s  Number 

5 

E 

TCA  (sec  from  A6A) 

6 

I 

L  =  Number  Data  Points 

station  file  is  terminated 

with  E.O.F. 

Solution  Parameters  and  Multisateilite  Covariance  (NCOV  Per  Tape) 

A.  Solution  Parameters 

Word  Number 

Format 

Value* 2  * 

1 

E 

Time  (A8-RSF)  in  sec 
(Relative  to  A6a) 

2 

E 

X  for  Satellite  1 

3 

E 

Y  for  Satellite  1 

4 

E 

Z  for  Satellite  1 

5 

E 

X  for  Satellite  1 

6 

E 

Y  for  Satellite  1 

7 

E 

Z  for  Satellite  1 

8 

E 

K,  —  Solar  Radiation  Parameter 
Satellite  1 

9 

E 

K2  -  Solar  Radiation  Parameter 
Satellite  1 

(2)Note:  Default  Value  is  9.9999E  +  9 
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Solution  Parameters  and  Multisatellite  Covariance  (NCOV  Per  Tape)  (Continued) 


Word  Number  Format 

10  E 

11  E 

12  E 

13  E 

Words  2-13  are  repeated  for  all  other 

12N  +  2  E 

12N  +  3  E 

12N  +  4  E 

12N  +  5  E 

12N  +  6  E 

12N  +  7  E 

•  • 

•  • 

•  • 

15N-  1  E 

15N  E 

15N+1  E 

15N  +  2  E 

15N  +  3  E 

15N+4  E 


Value 


K3  -  Solar  Radiation  Parameter 
Satellite  1 

Tj  -  Thrust  Parameter  Satellite  1 
T2  -  Thrust  Parameter  Satellite  1 
T3  -  Thrust  Parameter  Satellite  1 

satellites. 

At  —  Time  Bias  for  Satellite  N, 

At  -  Time  Drift  for  Satellite 
At  -  Aging  Rate  for  Satellite  N, 
At  -  Time  Bias  for  Satellite  N2 
At  -  Time  Drift  for  Satellite  N2 
At  -  Aging  Rate  for  Satellite  N2 


At  -  Time  Bias  for  Satellite  Nn 
At  -  Time  Drift  for  Satellite  Nn 
At  -  Aging  Rate  for  Satellite  Nn 
At  —  Time  Bias  for  Station  1 
At  -  Time  Drift  for  Station  1 
At  -  Aging  Rate  for  Station  1 


15N  +  3M-  1  E 

15N  +  3M  E 

15N  +  3M  +  1  E 

15N  +  3M  +  2  E 

15N  +  3M  +  3  E 

15N  +  3M  +  4  E 


At  Time  Bias  for  Station  M 
At  Time  Drift  for  Station  M 
At  Aging  Rate  for  Station  M 
a  priori  Sigma  on  At 
(Master  Station  m) 
a  priori  Sigma  on  At 
(Master  Station  m) 
a  priori  Sigma  on  At 
(Master  Station  m) 
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B.  Considered  Parameters 


1 

2 

3 

4 

5 

6 


E 

E 

E 

E 

E 

E 


X  -  Station 

Y  -  Station 
Z  —  Station 
X  -  Station 

Y  -  Station 
Z  -  Station 


t  (assumed) 
1  (assumed) 

1  (assumed) 

2  (assumed) 
2  (assumed) 
2  (assumed) 


3M  -  2  E 

3M-  1  E 

3M  E 

3M+  1  E 

3M  +  2  E 

3M  +  3  E 

3M  +4  E 

3M  +  5  E 


X  -  Station  M  (assumed) 

Y  -  Station  M  (assumed) 

Z  -  Station  M  (assumed) 

CR  -  Refraction  Connection 
Assumed  (=  0.0) 
a  priori  Sigma  (m)  on 

1,  4,  7,  etc. 

a  priori  sigma  (m)  on 

2,  5,  8,  etc. 

a  priori  Sigma  (m)  on 

3,  6,  9,  etc. 

a  priori  sigma  (%)  on 
3M  +  1 


C.  Multisatellite  Covariance  Upper  Triangular  Halt 
l  Row  Per  Record  Total  12N  +  3(N  +  M (Records 


Word  Number 


Format 


Definition 


1 

E 

Cov 

XI. 

XI 

( 1 .  1 )  Element  ' 

2 

E 

Cov 

XI, 

Y1 

(1,  2)  Element 

3 

E 

Cov 

XI, 

Z1 

(1.  3)  Element 

4 

E 

Cov 

XI, 

XI 

(1.4)  Element 

5 

E 

Cov 

XI, 

Y1 

(1,  5)  Element 

Dynamical 

6 

E 

Cov 

XI, 

Z1 

(1,  6)  Element 

Parameters 

7 

E 

Cov 

XI, 

K1 

(1.  7)  Element  , 

)  Satellite  1 

8 

E 

Cov 

XI, 

K2 

(1,  8)  Element  I 

9 

E 

Cov 

XI, 

K3 

(1,  9)  Element  ' 

i 

10 

E 

Cov 

XI, 

T1 

(1,  10)  Element 

1 

11 

E 

Cov 

XI, 

T2 

(1 ,  11)  Element 

| 

12 

E 

Cov 

XI, 

T3 

(1,  12)  Element 

f 

13 

E 

Cov 

XI, 

X2 

(etc.)  Element 

i 
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C.  Multisatellite  Covariance  Upper  Triangular  Halt  (Continued) 


Word  Number 

Format 

Definition 

• 

• 

• 

J  Dynamical 

• 

• 

• 

i  Parameters 

* 

• 

• 

'  Satellite  2 

24 

E 

Cov  XI,  T3  \ 

(N-  1)12  +  1 

E 

Cov  XI,  Xn 

(N-  1)12  +  2 

E 

Cov  XI,  Yn  ' 

,  Dynamical 

• 

• 

1 

Parameters 

• 

:  J 

Satellite  N 

(N)  12 

E 

Cov  XI,  T3 

E 

E 

Cov  XI,  TBj  ] 
Cov  XI,  TB,  1 
Cov  xi,  tbJ 
Cov  XI,  TB2] 

Time  Satellite 

E 

• 

• 

Cov  XI,  TB2  1 

Time  Satellite  2 

E 

Cov  XI,  tb2J 

• 

Cov  XI,  TBn1 

• 

I 

Time  Satellite  N 

E 

Cov  XI,  tbn 

E 

Cov  XI,  TB, 
Cov  XI,  TB,  1 
Cov  XI,  TB,  ] 
Cov  XI,  TBm 

E 

Time  Station  1 

• 

• 

• 

• 

■  Time  Station  M 

• 

• 

(Master  Station) 

E 

Cov  X1,TBm  , 

Row  2 

1 

E 

Cov  Y 1  YI 

2 

E 

Cov  Yl  ZI 

3 

E 

Cov  Yl  XI 

etc. 

etc. 

Row  3 

1 

E 

Cov  Zl  Zl 

2 

E 

Cov  Zl  XI 

3 

E 

Cov  Zl  Yl 

etc. 

etc. 
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REQUIRED  INPUTS  AND  SUGGESTED  DATA  CARDS 


Data  Quantities 


♦User  (1)  =-  3HAPL 

♦User  (2)  =  10H 

♦Version  (2)  =  10H 

♦Creator  (1)  =  10  HNSWC/DAHLG 

♦Creater  (2)  =  1  OH  REN 

♦Version  (2)  =  10H 

Inputs  (Not  intended  to  be  a  complete  or  limiting  list) 

IDATE  -  (Yr  -  1900)  *  1000  +  Day  Number 

♦Version  ( 1 )- alphanumeric  variable  for  output  on  TED  File 

♦NSAT  -  number  of  satellites  to  be  processed  (maximum  =  10) 

♦NSTAT  -  number  of  stations  for  which  data  are  present  (maximum  number  =  7) 
♦NCOV  —  number  of  times  covariance  matrix  is  to  be  output  on  TED  (maximum 
number  =10) 

♦MSAT(I)  -  satellite  number  for  each  satellite  being  processed 

♦MSTAT(I)  -  station  number  for  each  station  (Note:  The  last  station  in  the  array  is 
always  considered  the  master  station) 

♦TCOV(I)  -  array  containing  times  to  which  covariance  matrix  is  to  be  propagated  and 
output  on  TED  file  (1=1,  10;  program  should  be  written  so  maximum 
value  of  I  can  be  changed  without  causing  problems) 

♦NDYNP  -  number  of  dynamic  parameters  in  solution;  to  be  used  for  computing  array 
sizes  and  location  of  variables  in  pass  matrices 
♦POLARX,  POLARY  -  x,  y  components  of  polar  motion  at  epoch 
*A5A,  A5B,  A6A,  A6B  -  (Yr  -  1900)  *  1000  +  Day  +  SEC  *  10- 5  times  from  Run 

Specification  Form  (RSF)  parts  A5A,  A5B,  A6A,  and  A6B 
RSTART,  RSTOP  —  begin  and  end  times  for  range  data  used  in  solution 
RDSTART,  RDSTOP  -  begin  and  end  times  for  range  difference  used  in  solution 
KEY  -  flag  to  indicate  whether  range,  range  difference,  or  both  types  of  data  are  to 
be  processed 

SIGR,  SIGRD  -  sigmas  used  on  noise  for  range  and  range  difference  data,  respectively, 
in  data  generation  runs 

RATR,  RATRD  -  data  rate  for  range  and  range  difference  data  used  in  data  generation 
runs 

♦SIGMA(I)  —  I  =  1,  27;  Sigma  values  for  use  in  solution.  Suggested  array  ordering: 
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so 


SIGM(l)  sigma  for  Ist  element  of  satellite  state  vector  (x) 

SIGM(2)  sigma  for  2nd  element  of  satellite  state  vector  (y) 

SIGM(3)  —  sigma  for  3rd  element  of  satellite  state  vector  (z) 

SIGM(4)  -  sigma  for  4th  element  of  satellite  state  vector  (x) 

SIGM(5)  sigma  for  5,h  element  of  satellite  state  vector  (y) 

SIGM(6)  -  sigma  for  6th  element  of  satellite  state  vector  (z) 

SIGM(7)  -  sigma  for  x  component  of  radiation  pressure 

SIGM(8)  sigma  for  y  component  of  radiation  pressure 

S1GM(9)  sigma  for  z  component  of  radiation  pressure 

SIGM(IO)  sigma  for  x  component  of  thrust 

SlGM(ll)  sigma  for  y  xomponent  of  thrust 

SIGM(12)  sigma  for  z  component  of  thrust 

S1GM(13)  -  sigma  for  satellite  clock  bias 

SIGM(14)  -  sigma  for  satellite  clock  drift 

SIGM(1 5)  —  sigma  for  satellite  clock  aging 

SIGM(16)  —  sigma  for  station  clock  bias 

SIGM(  17)  -  sigma  for  station  clock  drift 

S1GM(18)  —  sigina  for  station  clock  aging 

S1GM(19)  -  sigma  for  master  station  bias 

S1GM(20)  —  sigma  for  master  station  drift 

S1GM(2 1 )-  sigma  for  master  station  aging 

S1GM(22)  -  sigma  for  x  component  of  station  position  (bias  parameter) 

S1GM(23)  -  sigma  for  y  component 
SIGM(24)  -  sigma  for  z  component 

SIGM(25)  -  sigma  for  x  component  of  station  position  (dynamic  parameters) 

SIGM(26)  -  sigma  for  y  component 

SIGM(27)  -  sigma  for  z  component 

SIGM(28)  -  sigma  for  refraction 

SIGM(29)  -  sigma  for  frequency  bias 

SIGM(30)  -  sigma  for  frequency  drift 

Most  variable  names  used  above  are  not  intended  to  restrict  the  programmer  in  his  own 
selection;  however,  it  would  be  convenient  if  those  variable  names  marked  with  *  were 
retained  as  they  frequently  occur  in  conversation  at  NSWC  and  in  communication  between 
NSWC  and  JHU/APL. 

Input  Files 

1.  Perturbed  Trajectory  Files  (1  for  each  satellite) 

2.  Pass  Matrix  Files  ( 1  for  each  data  class  of  each  satellite) 

3.  Polar  Motion  File  (see  Appendix  F  for  format) 

4.  ABCD  File  (see  Appendix  G  for  format) 
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appendix  c 

description  of  celest  files 
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PROPAGATED  TRAJECTORY  FILE  CEARTH-FIXEO  OR  INERTIAL) (TAPE51) 
THE  PASS  HATRIX  FILE  CONTAINS  THE  B-MATRIX  FOR  EACH  PASS 
OF  DA i A  INCLUDED  IN  THE  SPAN  OF  THE  T RA JECT ORY -  USED  FOR  FILTER¬ 
ING  THE  OATA.  ONLY  THE  ACCEPTED  PASSES  ARE  WRITTEN  ON  THE 
PASS  MATRIX  FILE.  RcCQRD  TYPE  i  IS  A  HEADER  RECORD  FOR  THE 
FILE  AND  APPEARS  ONLY  AS  THE  FIRST  RECORD  OF  EACH  FILE.  INFOR¬ 
MATION  I h  THIS  RECORO  IS  COMPARED  WITH  THE  SAME  TYPE  INFORMA¬ 
TION  FRO P  THE  TRAJECTORY  IN  THE  BSOLVR  SECTION  OF  CELEST.  ^IF 
THERE  iS  A  DIFFERENCE  IN  INFORMATION,  AN  ERROR  STOP  OCCURS. 
RECORD  TYPES  2  THRU  5  ARE  REPEATED  FOR  EACH  OATA  PASS. 

FORMAT  OF  PASS  MATRIX  FILE 

RECORD  TYPE  1 

WORD  SYMBOLIC 

NO.  TYPE  NAME  DESCRIPTION 


1  F  TRAT(i)  YEAR  OF  EPOCH  OF  TRAJECTORY 

2  F  THAT (2)  DAY  . 

3  A  I  TIME  TIME  CLOCK  VALUE  WHEN  TRAJECTORY  HAS  MADE 

4  I  10(1)  SATELLITE  NUMBER 

5  A  SAT (1)  SATELLITE  NAME 

6  A  SAT (2)  SATELLITE  NAME  (CONT) 

7  I  ID  < 21  SATELLITE  NUMBER  OF  2N0  SATELLITE 

8  A  SAT2(1)  SATELLITE  NAME  OF  2N0  SATELLITE 

9  A  SAT2 (2)  SATELLITE  NAME  (CONT) 
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53 


KECORO  TYPE  2 


WORD  SYMBOLIC 

NO.  TYPE  NAME  DESCRIPTION 


1  I  IPASS 


2  F  YEAR 

3  F  OAY 

4  A  ISTA 

5  I  NO 

6  I  ITYPE 

7  I  NO 

8  I  ITOLP 

9  I  NOF 

10  I  IPS 

11  0  PL 


12  F  W 


PASS  NO.  -  EACH  OATA  PASS  READ  FROM  THE 

OBS  FILE  IS  COUNTED.  IF  PASSES  ARE  REJECTED 
THE  PASS  NOS.  41 LL  NOT  3E  CONSECUTIVE. 

YEAR  OF  OBSERVATION  PASS 
DAY  OF  OBSERVATION  PASS 
STATION  NO.  OF  OBSERVING  STATION 
DATA  CLASS 
OATA  TYPE 

NUMBER  OF  OBSERVATIONS  IN  PASS 
TOTAL  NUMBER  OF  PARAMETERS  IN  B-MATRIX, 
INCLUDING  BIASES 
NO.  OF  ACCEPTEO  POINTS  IN  PASS 
PASS  STATUS  =0. GOOD, NOT  = 0 ,bAD ( ALWAYS  O.NOW) 
PARAMETER  LABEL.  THIS  WORD  INDICATES  WHICH 
PARAMETERS  MAKE  UP  THE  MATRIX.  A  BINARY  BIT 
IS  ALLOWED  FOR  EACH  PARAMETER  LEFT  AOJUSTEO 
IN  THE  FOLLOWING  OROER.  X,  Y,  Z,  XDOT,  Y  DOT , 
ZDOT,  OR  (  A,  ES,  EC,  I,  AMO,  OMEGA),  CO,  FX, 
FY, FZ,KR,XS,YS,ZS, REF, F, FOOT 
VARIENCE 
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kF  Co k 0  TYPE  2 


(  cour.) 


WGk  0 

SYMBOLIC 

NO. 

TYPt 

NAME 

DESCRIPTION 

13 

F 

ELEV 

ELEVATION  AT  TCA  OF  PASS 

( DEG ) 

14 

F 

PHIO 

STATION  LONGITUDE  lUEGi 

15 

*  F  ’ 

FNS 

filtered  noise 

16 

F 

FSS 

SATELLITE  FREUUEMCY 

17 

F 

x  :■ 

EARTH-FIXED  X  COMPONENT 

OF 

station 

POSITION 

18 

F 

YO 

EARTH-FIXED  Y  COMPONENT 

OF 

STATION 

POSITION 

19 

F 

z , 

earth-fixeo  z  component 

STATION 

POSH  ION 
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RECORD  TV  Ft  3 


WORD  SYMBOLIC 

NO.  '  VPE  NAME  DESCRIPTION 


1  F  dM  AT  LONER  T  RI ANGULAR  PART  OF  B-EATRIX  FOR  PASS 

2  F 


*N  F 

N+l  F  EtfEC(l)  E-VECTOR  FOR  PASS 
N  +  2  F 


•  • 

•  • 

♦♦  M  F 


EwEC  (ITOLPI 


♦WHERE  N  =(ITOLP*(ITOLP+HJ/2 
♦♦WHERE  M  •  N+ITOLP 
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RECORD  TYPE  4 


MORO  SYMBOLIC 

NO.  TYPE  NANE  DESCRIPTION 


1  I  NOdlAS  NO.  OF  9IASES 

2  I  IHRMIN  RISE  TINE  OF  PASS  CPACKEO  AS  HR  ANO  MIN) 

*  3  f  IOUR  '  OURAT IOfl  OF  PASS  (MlN) 

4  I  I  AC  AUTOCORRELATION 

5  I  I VERS  INDICATES  FROM  HHI CH  CELEST  VERSION  FILE  MADE 

6  F  TRISE  RISE  TINE  OF  PASS  CSEC.  FROM  MIDNIGHT* 

7  F  TSET  SET  TIME  OF  PASS  CSEC.  FROM  MIDNIGHT! 


RECGRO  TYPE  5 


MOkO  SYMBOLIC 

NO.  TYPE  NAME  DESCRIPTION 


1  F  TCA  TIME  OF  CLOSEST  APPROACH  3F  PASSISEC.  FROM 

MIDNIGHT) 

2  F  RM AT (1 » 1 )  ROTATION  MATRIX  TO  ROTATE  FFOM  INERTIAL 

3  F  RMATC1.2)  TO  EARTH-FIXED  AT  TCA. 


•  .  • 

10  F  RMAT  <3*31 


RECORD  TYPES  2  THRO  5  ARE  REPEATED  FOR  EACH  DATA  PASS 
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INERTIAL  PEkTUR3tU  TRAJECTORY  FKE  (TA=>E9I 

THE  INERTIAL  PEETURbED  TRAJECTORY  FILE  CONSISTS  OF  TNO  RECORD 
TYPES  -  RE CORO  TYPE  1  BEING  A  HEADER  RECORD  APPEARING  ONLY  ONCE, 

ANU  RECORD  TYPE  2  WHICH  CONTAINS  THE  POSITION  OF  THE  SATELLITE 
ANO  PARTIALS  OF  POSITION  WITH  RESPECT  TO  EACH  PARAMETER. 

ONLY  1  HE  CANONICAL  PARTIALS  WRT  ORAG  ANO  THRUST  APPEAR  ON  THE  TRAJ. 
THE  PARAMETER  OROER  IS  AS  FOLLOWS.  THE  SIX  0R3IT  PARAHETES • ARE  FIRST, 
EITHER  COORDINATES  OR  ELEMENTS,  FOLLOWED  BY  DRAG,  THREE  THRUST,  ANO 
RAOIATION  PRESSURE.  DRAG,  THRUST,  ANO  RADIATION  ARE  OPTIONAL.  A 
SYMBOLIC  LIST  ANO  OROER  OF  PARAMETERS  IS  AS  FOLLOWS.  X,  Y,  Z,  XOOT, 
YOOT,  ZOOT,  OR  (  A,  ES,  EC,  I*  AMO,  OMEGA),  CO,  FX,  FY,  FZ,  KR 
RECOku  TYPE  2  IS  REFEATEO  FOR  EACH  TINE  LINE  0*  THE  TRAJECTORY. 


FORMAT  OF  INERTIAL  PERTURBED  TRAJECTORY  FILE 


RECORD  1 YPE  i 


WORO 

NO. 

1 

2 

3 

4 

5 

6 
7 
6 

% 

10 

11 

12 

13 

14 


SY  M JOLIC 

TYPE  NAME  OESCRIFTION 


1 

F 

F 

F 


NB  RECORO  NO.  =  1 

TRAT (1)  YEAR  OF  EPOCH  OF  TRAJECTORY 
TAAT<2)  OAY 

TRAT  <31  SEC  “ 


F 

F 

F 

F 

1 

I 

1 

1 

A 

X 


TRAT <41  TVE-TO,  TIME  UF  V £ RNAL  EQUINOX  MINUS  EPOCH 
OF  THE  TRAJECTORY 

TRAT  <51  INTERVAL  AT  WHICH  TRAJ  IS  WRITTEN  (SEC) 
TRAT (6)  LAST  TIME  ON  TRAJECTORY (S£C  FROM  EPOCH) 
TRAT <7 )  INTEGRATION  INTERVAL  USED  WHEN  CREATING 

THE  TRAJECTORY  (SEC) 

IFLOH  KINO  OF  TRAJECTORY  =  6, 7, 4,0*  5, 

10(1)  SATELLITE  NUMBER 

KCG(IO)  INOICATES  IF  TRAJECTORY  WAS  MAOE  USING 

EPOCH  OF  OATE= 1,  OR  EPOCH  OF  1950=  0 
10(3)  IMPROVEMENT  CYCLE  NUMBER 

IOATE  TIME  CLOCK  VALUE  WHEN  TRAJECTORY  WAS  HADE 

NPAR  NUMBER  OF  PARAMETERS 
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RECORD  NO.  1  <  CONT.) 

WOr 0  SYMBOLIC 

NO.  TYPE  NAME  DESCRIPTION 


15  0  LAdELP  PARAMETER  LABEL  WORO.  THIS  WORD  INDICATES 

WHICH  PARAMETERS  WERE  USED  WHEN  THE  FILE  WAS 

•  •  •  •  4  •  9  •  •  r  .  , 

MAOE.  A  BINARY  BIT  IS  ALLOWED  FOR  EACH  PARA¬ 
METER  INDICATED  A o OWE.  LErT  ADJUSTED.  FOR 
EXAMPLE*  X  OR  A  IS  INDICATEC  IN  THE  LEFT  MOST 
BIT  OF  THE  WORO,  Y  OR  ES  IS  IN  THE  NEXT.  ETC. 
IF  THE  BIT  IS  ON  (  =  1).  THEN  THE  PARAMETER  MAS 
USED.  IF  THE  BIT  IS  OFF  (=0),  THE  PARAMETER 
WAS  NOT  USED. 


16 

I 

NOPA 

NO. 

OF  ORAGS  USED  WHEN  MAKING 

TRAJECTORY 

17 

I 

NT  PA 

NO. 

OF  THRUST  USED  WHEN  MAKING  THE 

TRAJ 

18 

F 

QT  I M  ( 1 ) 

END 

TIME  OF  1ST  DRAG  SEGMENT 

(SEC 

FROM  EPOCH) 

19 

F 

0TIM(2I 

•• 

-  **  2ND 

•I 

M  — 

• 

F 

• 

«• 

H  ••  ••  •• 

M 

•«  M 

• 

F 

• 

«• 

M  ••  •• 

•• 

N  •• 

• 

F 

• 

•• 

N  ••  II  M 

•• 

N  •• 

N 

F 

DT I  MS  CNDPA ) 

END 

TIME  OF  LAST  DRAG  SEGMENT 

WHERE  N=(>GPAU7  OR  N=18*  IF  NOPA=3 


C-ll 


u 


RECORD  NC» 

1  (  CONT 

.) 

MOKO 

SYMBOLIC 

NO* 

TYPE 

NAMt 

DESCRIPTION 

— 

----- 

— 

N*1 

F 

TT1NS  (1*1) 

START 

TIME  OF  1ST 

THRUST  SEGMENT 

N+2 

F 

TTIMS(2»i) 

ENO 

H  M  •• 

••  •• 

• 

F 

• 

•  •  •  • 

• 

* 

F 

• 

• 

F 

• 

• 

F 

• 

♦H-l 

F 

TTIHS(l.L) 

iT  ART 

TINE  OF  LAST 

THRUST  SEGMENT 

M 

F 

TTIMS(2.L) 

END 

M  — 

M  +  l 

F 

01(1) 

X  ) 

M+2 

F 

01(2) 

Y  I  • 

INERTIAL  COMPONENTS  OF  SATELLITE 

M*3 

F 

01(3) 

Z  ) 

POSITION  AT  EPOCH  (KM) 

M+L 

F 

01(4) 

XOOT 

) 

M  +  5 

F 

01(5) 

YOOT 

f-  INERTIAL 

COMPONENTS  OF  SATElITE 

M+6 

F 

01(6) 

ZOOT 

)  VELOCITY  AT  EPOCH  (KM/SEC) 

M  +  7 

F 

01(7) 

C01  - 

COEFFICIENT 

OF  DRAG  FOR  1ST  DRAG 

SEG 

M*cj 

F 

01(6) 

CD  2  - 

• 

,  .  .  2ND  . 

• 

M*9 

F 

01(9) 

C03  - 

• 

,  •  «  3R0  • 

• 

H+ir 

F 

01(10) 

CD  4 

• 

,  .  .  4TH  . 

• 

M+ll 

F 

01(11) 

CO  5 

• 

.  .  .  5TH  . 

• 

M+12 

F 

01 (12) 

CD  6 

• 

.  .  .  6TH  . 

• 

M  +  13 

F 

01(131 

CO  7 

• 

• 

• 

• 

7TH 

M+14 

F 

01 (141 

CO  8 

• 

• 

• 

• 

8TH 

M*15 

F 

01(151 

CD  9 

• 

• 

• 

9TH 

M+16 

F 

01 (16) 

C010 

• 

• 

• 

• 

10TH 

M+17 

F 

01 (17) 

CD  11 

• 

• 

• 

• 

11TH 

M*18 

F 

01 (18) 

CO  12 

• 

• 

• 

• 

12TH 

H+19 

F 

01 (19) 

CO  13 

• 

• 

• 

• 

13TH 

M*20 

F 

01  (20 ) 

C014 

• 

• 

• 

• 

14TH 

M+21 

F 

01(21) 

CO  15 

• 

• 

• 

• 

15TH 

M+22 

F 

01(22) 

CO  16 

• 

• 

% 

• 

16TH 

M*23 

F 

01(231 

CO  17 

• 

# 

• 

• 

17TH 

M+24 

F 

01(24) 

C018 

• 

• 

• 

• 

18TH 

M+25 

F 

01(25) 

CO  19 

• 

• 

• 

« 

19TH 

M-*26 

F 

01(26) 

CO  20 

• 

• 

• 

• 

20TH 

WHERE 

M= 

N+2MNTPA) 

UR  M=N*2, 

IF  NT PA= 

C-13 


RECORD  NC.  1 


(  CONT.I 


woo 

SYMBOLIC 

•  1 

O  1 

Z  1 

TYPE 

NAME 

DESCRIPTION 

M+27 

F 

01(27) 

AX  1 

) 

M+2a 

F 

01 (261 

AY  1 

(-COMPONENTS  OF  1ST  IMPOST (KM/SEC**2) 

M+29 

F 

01(29) 

AZ1 

) 

M+3l 

F 

01(30) 

AX  2 

) 

M+31 

F 

Oi (31) 

AY  2 

1 - COMPONE  NTS  OF  2ND  THPUST (KM/SEC**2) 

M+32 

F 

01(32) 

AZ2 

) 

M+33 

F 

01 (33) 

AX  3 

) 

M+34 

F 

01(34) 

AY  3 

(-COMPONENTS  OF  3R0  THRUST (KM/SEC**2) 

M+35 

F 

01 (35) 

AZ3 

) 

M+36 

F 

01(36) 

AX  4 

) 

M+37 

F 

01  (37) 

AY  4 

(-COMPONENTS  OF  4TM  T HPUST ( KM/SEC** 2) 

M+38 

F 

01(38) 

AZ4 

M+39 

F 

01 (39) 

KR 

-  RADIATION  PRESSURE  COEFFICIENT 

M+40 

F 

01(40) 

KR2 

-2ND  RADIATION  PPESSURE  COEFFICIENT 

M+41 

F 

01 (41) 

KR3 

-3RD  RAOIATION  PRESSURE  COEFFICIENT 

NGPA 

NO. 

OF  GRAVITY  PARAMETERS 

10G0 

GRAVITY  ID  ARRAY  (32) 

OLPGP(l) » 

II 

♦NGPA  DELTA  P  GRAVITY 

NRPA 

NO 

OF  RAOIATION  PARAMETERS  BEING  IMPROVEO 
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RECORD  1 Y  Ft  2 


wOt  0 

SYMBOLIC 

NO.  TYPE 

NAME 

DESCRIPTION 

1 

I 

N8 

CONSE 

CUTIVE  RECORD  NO. 

=  2.3, 

4*.  .  N 

2 

F 

TI 

TI 

-TO 

(SECONOS  FROM  EPOCHI 

3 

F 

TRA(2) 

X 

) 

4 

F 

TRA<3) 

Y 

!  - 

INERTIAL  COMPONENTS  OF 

SATELLIT 

£ 

5 

F 

TRA ( 4) 

Z 

I 

POSITION  AT  TIME 

(TI) 

(KM) 

6 

F 

.TRA^l 

PARTIAL  X 

AT  (TI)  WRT 
• 

PARAMETER  (1) 

AT 

(TO) 

7 

F 

TRA(6) 

• 

Y 

t  »  • 

• 

(1) 

• 

• 

8 

F 

TRA  (7 ) 

• 

Z 

•  •  • 

• 

(1) 

• 

• 

*3 

F 

TkA(c) 

• 

X 

•  •  • 

• 

(2) 

• 

• 

10 

F 

TkA(9) 

• 

Y 

•  •  • 

• 

(2) 

• 

• 

11 

F 

TkA(lC) 

t 

Z 

•  •  • 

• 

(2) 

• 

• 

« 

• 

• 

• 

• 

•  •  • 

• 

• 

• 

• 

n 

• 

• 

• 

• 

•  •  • 

• 

• 

• 

• 

• 

• 

• 

• 

• 

•  •  • 

• 

• 

• 

• 

N-2 

F 

TRA(M-2I 

• 

X 

•  •  • 

• 

NPAR 

• 

• 

N-l 

F 

TxA(M-l) 

• 

Y 

•  •  • 

• 

NPAR 

• 

• 

N 

F 

TRA (M) 

• 

Z 

•  •  • 

• 

NPAR 

• 

• 

WHERE  NEAR  =  NO.  OF  PAkAMETERS 

=  6  ORBIT  +  MIN0  (l,NOPA)+3*MlNb  (i.NTPA) +MINO (l.NFPA) 
nOPA-NO.  OF  ORAG  SEGMENTS 
N; pa-no.  of  thrust  segments 

NrPA=NO.  OF  RADIATION  PRESSURE  PARAMETERS  <0  OR  -II 
N-  5  +3*NP AR 
M-  4+3*NFAR 
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0IM6N0STIC  INFORMATION  FILE  (TmPEISI 

THE  DIAGNOSTIC  INFORMATION  FILE  IS  CREATED  BY  THE  COWAR 
SECTION  OF  CElEST  TO  BE  USED  BY  THE  GRAPHICS  OR  POST  ANALYSIS 
SECTION  FOR  CELEST.  IT  CONTAINS  STATISTICAL  (3  AT  A  WHICH  IS 
USED  AS  A  MEASURE  OF  HON  GOOD  THE  FIT  OVER  A  PARTICULAR  SPAN  IS. 
IT  CONSISTS  OF  FOUR  RECORD  TYPES.  RECQRO  TYPES  1,2.  ANO  A  APPEAR 
ONLY  ONCE  FOR  A  REV.  oUT  ARE  REPEATED  IF  MORE  THAN  ONE  REV  IS 
PkOCESSEO.  RECORD  TYPE  3  IS  REPEATED  FOR  EACH  OUTPUT  DELTA  T 
OUR I  No  THE  SPAN  OF  THE  REV. 

FORMA!  OF  DIAGNOSTIC  INFORMATION  FILE 


RECORD  TYPE  1 

WOi  0  SYMBOLIC 

NO.  TYPE  NAME 


1  F  TO (II 

2  F  T042I 

3  F  TO  <31 

4  I  10(1) 

5  I  CS 

6  I  NREV 

7  F  TK AT ( II 

8  F  T  RAT ( 21 

9  F  TRAT (31 

1C  F  TDlT 

11  I  MT 

RECURO  TYPE  2 

NO.  0  SYMBOLIC 

NO.  TYPE  NAME  DESCRIPTION 


1  F  IlT  LAST  TIME  FOR  WHICH  RECORD  TYPE  3  IS 

WRITTEN  (SEC.  FROM  REV  EPQCHI 


DESCRIPTION 


YEAR  OF  EPOCH  OF  REV 

OA  Y  .  .  .  . 

SEC.  .... 

SATELLITE  NUMBER 

SCALE  FACTOR  FOR  COVARIAMDE 

KEV  NUMBER  OF  SATELLITE 

YEAR  OF  EPOCH  OF  F ERT  TRAJ  (LONG  ARCI 

DAY  ...  .  .  .  * 

SEC.  ...  .  .  a  • 

OUTPUT  OELTA  T  FOR  EARTH-FIXED  ‘VALUES 
MULTIPLE  OF  TDLT  TO  OUTPUT  INERTIAL  VALUES 
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RECuKO  TYRE  3 


HGnD  SYMBOLIC 

NO.  YPE  NAME  DESCRIPTION 


1  (r  JIM  TIME  IN  SEC.!  F  RDM  REV  EPOCH 

2  F  oR(itl)  )  CONFI3ENCE  IN  THE  TRAJECTORY 

3  F  8K(2t2)  ! -  RAOIAL*  TANGENTiALt  ANQ  OUT 

4  F  8R (3. 3)  )  DIRECTION  AT  TIME  (TIM) 

RECORD  TYPE  3  IS  REPEATED  FOR  EACH  OUTPUT  TIME  IN 


IN  THE 
OF  PLANE  (R 

REV  SPAN. 


THE 

LAST 

TYPE  3 

RECORD  CONTAINS  FOUR  WORDS  OF  ZERO  VALUE. 

RECO 

-  C  1  Y 

PE  4 

WUkO 

SYMBOLIC 

NO. 

TYPE 

NAME 

DESCRIPTION 

1 

F 

Si  G2R 

AVERAGE  SIGMA  RADIAL  FOR  ENTIRE  REV 

2 

F 

SIG2V 

AVERAGE  SIGMA  TANGENTIAL  FOR  ENTIRE  REV 

3 

F 

S2RXV 

AVERAGE  SIGMA  <RXV  )  FOR  ENTIRE  REV 

i 


appendix  d 

SUBROUTINE  NTRP8 
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SUBROUTINE  NTRP8 


Common/NTRP/IPRNT,  ITRAJ,  IR,  PCTIO,  1TRA1,  1TRA2,  TRA(8,133),  TRAK267) 
TRAT(7),  CTIO 

TRAT(  1 )  =  year  of  Epoch  "j 

TRAT(2)  =  day  of  Epoch  I 

TRAT(3)  =  seconds  of  Epoch  ?  from  header  record  of  trajectory  file 
TRAT(4)  =  tvc  -  t0  (GHA)  I 

TRAT(5)  =  Atw  ; 


IPRNT  -  output  file,  so  error  messages  can  be  written 
ITRAJ  -  intput  trajectory  (perturbed  trajectory)  file 

TRA  -  array  containing  records  2-9  of  trajectory  file;  these  must  be  read  in  before  lsl  call 
to  NTRP8 

TRAI  -  array  containing  the  interpolated  values 

TRAI(l)-x  TRAI(7)  —  x 

TRAI(2)  -  y  TRAI(8)  -  y 

TRAI(3)  -  z  TRAI(9)  -  z 

TRAI(4)  -  x 
TRAI(5)  -  y 
TRAI(6)  -  z 

NPAR  -  number  of  parameters  (from  header  record) 

CTIO  —  time  in  seconds  for  which  to  interpolate  (from  A5A) 

PCTIO  -  previous  time  for  which  interpolation  was  done;  initialize  to  -9999.  in  calling 
program 

ITRA1  -  number  of  words  in  last  record  read  from  trajectory  file  (NPAR  *3 +4) 

ITRA2  -  block  number  of  last  record  read  from  trajectory  file 

IER  -  error  return  from  NTRP8 
=  0  interpolation  was  done 

=  1  time  to  interpolated  occurred  too  early  to  interpolate  (before  fourth  time  on 
trajectory) 
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=.  2  time  to  interpolate  occurred  too  late  to  interpolate  (not  enough  data  on 
trajectory) 

1R  -  indicates  what  to  interpolate  for 
=  1,  interpolate  for  r,  r 
=  2,  interpolate  for  r,  r,  r 
=  3,  interpolate  for  i px 
=  4,  interpolate  for  i//x,  ^ 

NOTE: 

1.  Subroutine  must  be  called  twice  (JR  =  1  and  1R  =  4)  when  position,  velocity, 
and  partials  are  required. 

2.  Subroutine  must  be  reinitialized  for  each  satellite  (trajectory  files  must,  of  course, 
be  rewound  before  reading  header  and  1st  records  prior  to  calling  NTRP8). 

3.  The  order  of  the  partials  in  the  TRA1  array  for  ipx  or  ip-  is 

partials  with  respect  to  position 
partials  with  respect  to  velocity 

partials  with  respect  to  drag  (not  used  by  the  SATRACK  program) 
partials  with  respect  to  three  components  of  thrust  (if  present  in  the  solution) 
partials  with  respect  to  three  components  of  radiation  pressure  (if  present  in  the 
solution) 

partials  with  respect  to  gravity  parameters  (if  present  in  the  solution) 


i 
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POLAR  POSITION  FILE 


E-l 


POLAR  POSITION  FILE 


REF;  “Polar  Position  File  System,”  R.  G.  Barker,  January  1975 

There  are  nine  records,  also  referred  to  as  subfiles,  indexed  by  integer  labels: 

1.  IPMS  -  preliminary  pole  position 

2.  IPMS  —  final  pole  position 

3.  BIH  -  preliminary  pole  position 

4.  BIH  -  final  pole  position 

5.  BIH  UT1-UTC 

6.  USNO  UT1-UTC 

7.  Bi-daily  solutions  from  the  CELEST  program 

8.  Predicted  pole  position 

9.  Predicted  UT1-UTC 

Data  in  each  subfile  are  in  order  by  day  of  the  year.  Space  is  reserved  for  each  day, 
whether  or  not  data  are  available  for  that  day.  The  first  word  is  always  year  number.  Only 
one  additional  word  per  day  is  needed  in  each  of  the  subfiles  for  UT1-UTC  (numbers  5.  6, 
and  9).  Two  words  are  needed  in  each  of  the  other  subfiles  for  the  x  and  y  coordinates 
of  polar  position.  Moreover,  two  words  are  allocated  for  the  standard  deviation  of  these 
coordinates  in  subfiles  4  and  7.  The  later  subfile,  namely  the  Bi-daily  solutions,  also 
provides  for  a  word  of  alphanumeric  characters  identifying  the  satellite  or  source  of  data. 
It  differs  from  the  other  subfiles  in  another  way:  three  complete  sets  of  data  may  be 
stored  each  day,  corresponding  to  different  satellites.  Thus,  18  words  are  set  aside  for 
every  day  in  the  subfile.  The  following  figure  indicates  the  various  subfile  structures.  All 
numerical  data  are  in  floating  point  format. 


Subfiles 


5,  6,  9 

year 

UT1-UTC 

1,  2,  3,  8 

year 

X 

y 

4 

year 

X 

y 

°x 

°y 

7 

year 

X 

y 

®x 

°y 

satellite 
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APPENDIX  F 

TRACKED  EPHEMERIS  DATA  TAPE  BINARY  FORMAT 
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TRACKED  EPHEMERIS  DATA  TAPE 
BINARY  FORMAT 


Revised  1/19/77 


Tape  Header  (1  Per  Tape) 

Word  Number  Format  Definition 


1,2  A 

3,4  A 

5,6  A 

7  I 

8  I 

8+1  1 

8  +  2  1 

• 

• 

8  +  N  I 

8  +  N+  1  I 

8  +  N  +  2,  8  +  N  +  3  A 

8  +  N+4  1 

8  +  N  +  5  E 

8+N+6  E 

8  +  N  +  7  E 

8  +  N  +  8  E 


APL  -  Tape  User 
NSWC/Dahlgren  -  Tape  Generator 
MODX  -  OD/CT  Version  X 
(Yr  -  1900)  *  1000  +  Day  Number 
Run  Specified  (A1  -  RSF) 

N  -  Number  Satellites 
N1  -  Satellite  Number 
N2  -  Satellite  Number 


NN  -  Satellite  Number 
M  -  Number  Stations 
RW  or  FW  -  Real  or  Fitted 
Trajectories 

NCOV  -  Number  of  Covariance 
Matrices  Written  Out 
(Yr  -  1900)  *  1000  +  Day  Number 
+  Sec/ 10s  (A5a-RSF) 

(Yr  -  1900)  *  1000  +  Day  Number 
+  Sec/105  (A5b-RSF) 

(Yr  -  1900)  *  1000  +  Day  Number 
+  Sec/105  (A6a-RSF) 

(Yr  -  1900)  *  1000  + Day  Number 
+  Sec/105  (A6b-RSF) 
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u. 


Ephemeris  Files  (1  Per  Satellite) 
A.  Header  Record* 1  * 


Word  Number 

Format 

Definition 

1 

I 

(Yr-  1900)  *  1000  +  Day  of 

Ephemeris  Record 

2 

E 

t  -  Seconds  UTC  After  Midnite 
(A6a-RSF) 

3 

I 

Satellite  Number 

4 

E 

X  Satellite  Number  (A6a-RSF) 

5 

E 

Y  Satellite  Number 

6 

E 

Z  Satellite  Number 

7 

E 

X  Satellite  Number 

8 

E 

Y  Satellite  Number 

9 

E 

Z  Satellite  Number 

10 

E 

True  Right  Ascension  of  Greenwich 
at  Epoch  A6A-PSF  (Radians) 

11 

E 

X  -  Complete  Polar  Motion  at  Epoch 

12 

E 

Y  -  Complete  Polar  Motion  at  Epoch 

13 

E 

(1,  1)\ 

14 

E 

0,  2) 

| 

15 

E 

(1,  3)  1 

[  Elements  of  Transformation 

16 

E 

(2,  1) 

Matrix  From  Mean  1950.0 

17 

E 

(2,  2) 

)  Inertial  to  Earth-Fixed 

18 

E 

(2,  3) 

Geocentric  Equatorial  at 

19 

E 

(3,  1) 

\  Solution  Epoch  (A6a-RSF) 

20 

E 

(3,  2) 

21 

E 

(3,  3)  , 

1 

\ 


Note:  Header  record  written  for  time  A6a  for  all  satellites;  for  satellite  number  1,  it  is  repeated  NCOV  times  with  the  right 
ascension  of  Greenwich  and  the  transformation  matrix  updated  to  correspond  to  the  times  that  the  covariance  matrix  will  be  output 
in  part  IV. 
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B.  Ephemeris  Data  Records  ( 1  Per  Time  Line ) 


Word  Number 

Format 

Definition 

1 

E 

Time  (sec)  Relative  to 

Epoch  (A6A) 

2 

E 

X  -  Coordinate  (km )  in 

Mean  1950.0 

3 

E 

Y  -  Coordinate  (km)  in 

Mean  1950.0 

4 

E 

Z  -  Coordinate  (km)  in 

Mean  1950.0 

5 

E 

X  -  Coordinate  (km)  in 

Mean  1950.0 

6 

E 

Y  —  Coordinate  (km)  in 

Mean  1950.0 

7 

E 

Z  -  Coordinate  (km)  in 

Mean  1950.0 

8 

E 

At  -  Satellite  Clock  Error 

epheremis  file  terminates  with  E.O.F. 

Station  Files  ( 1  Per  Station) 

A.  Header  Record 

Word  Number 

Format 

Definition 

1 

I 

Station  Number 

2 

E 

0  -  Latitude  (deg)  in 

WGS-72  System 

3 

E 

X  -  Longitude  (deg)  in 

WGS-72  System 

4 

E 

h  -  Height  (m)  in  WGS-72  System 

5 

I 

Data  Class 

1  =  Range 

2  =  Range  Difference 

3  =  Both 

6 

E 

Mean  on  Range  Data  (=  0.0) 

7 

E 

Sigma  of  Noise  on  Range  Data 
in  Data  Generation  R 

8 

E 

Range  Data  Rate  -  Data 

Generation  Run  (sec) 
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01.  Station  Files  (1  Per  Station)  (Continued) 


Word  Number 

Format 

Definition 

9 

E 

Mean  on  Range  Difference  Data 
(=  0.0) 

10 

E 

Sigma  of  Noise  on  Range  Difference 
Data  in  Data  Generation  R 

11 

E 

Range  Difference  Data  Rate  - 
Data  Generation  Run 

12 

I 

K  -  Total  Passes  From  Station 

B.  Data  Records  ( 1  Per  Pass) 

1 

1 

Pass  Number 

2 

E 

Start  Time  (sec  from  A6a) 

3 

E 

End  Time  (sec  from  A6a) 

4 

I 

Tracked  Satellite’s  Number 

5 

E 

TCA  (sec  from  A6A) 

6 

I 

L  =  Number  Data  Points 

station  file  is  terminated  with  E.O.F. 

Solution  Parameters  and  Multisatellite  Covariance  (NCOV  Per  Tape) 

A.  Solution  Parameters 

Word  Number 

Format 

Value* 2  * 

1 

E 

Time  (A8-RSF)  in  sec 
(Relative  to  A6a) 

2 

E 

X  for  Satellite  1 

3 

E 

Y  for  Satellite  1 

4 

E 

Z  for  Satellite  1 

5 

E 

X  for  Satellite  1 

6 

E 

Y  for  Satellite  1 

7 

E 

Z  for  Satellite  1 

8 

E 

K,  -  Solar  Radiation  Parameter 
Satellite  1 

9 

E 

K2  -  Solar  Radiation  Parameter 
Satellite  1 

< 2)Note:  Default  Value  U  9.9999E  +  9 
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Solution  Parameters  and  Multisatellite  Covariance  (NCOV  Per  Tape)  (Continued) 


Word  Number  Format 

10  E 

11  E 

12  E 

13  E 

Words  2-13  are  repeated  for  all  other 

12N  +  2  E 

12N  +  3  E 

12N  +  4  E 

12N  +  5  E 

12N  +  6  E 

12N  +  7  E 

•  • 

•  • 

•  • 

15N-  1  E 

15N  E 

15N+1  E 

15N  +  2  E 

15N  +  3  E 

15N  +  4  E 


Value 


K3  -  Solar  Radiation  Parameter 
Satellite  1 

T,  -  Thrust  Parameter  Satellite  1 
T2  -  Thrust  Parameter  Satellite  1 
T3  -  Thrust  Parameter  Satellite  1 

satellites. 

At  -  Time  Bias  for  Satellite  N, 

At  -  Time  Drift  for  Satellite  N, 
At  -  Aging  Rate  for  Satellite  Nj 
At  -  Time  Bias  for  Satellite  N2 
At  -  Time  Drift  for  Satellite  N2 
At  -  Aging  Rate  for  Satellite  N2 


At  -  Time  Bias  for  Satellite  Nn 
At  -  Time  Drift  for  Satellite  Nn 
At  -  Aging  Rate  for  Satellite  Nn 
At  -  Time  Bias  for  Station  1 
At  -  Time  Drift  for  Station  1 
At  -  Aging  Rate  for  Station  1 


15N  +  3M-  1  E 

15N  +  3M  E 

15N  +  3M+1  E 

15N  +  3M  +  2  E 

15N  +  3M  +  3  E 

15N  +  3M  +  4  E 


At  Time  Bias  for  Station  M 
At  Time  Drift  for  Station  M 
At  Aging  Rate  for  Station  M 
a  priori  Sigma  on  At 
(Master  Station  m) 
a  priori  Sigma  on  At 
(Master  Station  m) 
a  priori  Sigma  on  At 
(Master  Station  m) 


B.  Considered  Parameter 


1 

2 

3 

4 

5 

6 


E 

E 

E 

E 

E 

E 


X  -  Station 

Y  -  Station 
Z  -  Station 
X  —  Station 

Y  -  Station 
Z  —  Station 


1  (assumed) 
1  (assumed) 

1  (assumed) 

2  (assumed) 
2  (assumed) 
2  (assumed) 


3M-  2  E 

3M-  1  E 

3M  E 

3M  +  1  E 

3M  +  2  E 

3M  +  3  E 

3M  +  4  E 

3M  +  5  E 


X  -  Station  M  (assumed) 

Y  -  Station  M  (assumed) 

Z  -  Station  M  (assumed) 

CR  —  Refraction  Connection 
Assumed  (=  0.0) 
a  priori  Sigma  (m)  on 

1,  4,  7,  etc. 

a  priori  sigma  (m)  on 

2,  5,  8,  etc. 

a  priori  Sigma  (m)  on 

3,  6,  9,  etc. 

a  priori  sigma  (%)  on 
3M+  1 


C.  Multisatellite  Covariance  Upper  Triangular  Halt 


1  Row  Per  Record  Total  12N  +  3(N  +  M)Records 


Word  Number  Format 


Definition 


1 

E 

Cov 

XI, 

XI 

0, 

1)  Element 

2 

E 

Cov 

XI, 

Y1 

(1, 

2)  Element 

3 

E 

Cov 

XI, 

Z1 

(1, 

3)  Element 

4 

E 

Cov 

XI, 

XI 

(1, 

4)  Element 

5 

E 

Cov 

XI, 

Y1 

(1, 

5)  Element 

6 

E 

Cov 

XI, 

Z1 

(1, 

6)  Element 

7 

E 

Cov 

XI, 

K 1 

(1, 

7)  Element 

8 

E 

Cov 

XI, 

K2 

(1, 

8)  Element 

9 

E 

Cov 

XI, 

K3 

(1, 

9)  Element 

10 

E 

Cov 

XI, 

T1 

(1, 

10)  Element 

11 

E 

Cov 

XI, 

T2 

(1, 

1 1 )  Element 

12 

E 

Cov 

XI, 

T3 

(1, 

12)  Element 

13 

E 

Cov 

XI, 

X2 

(etc.)  Element 

Dynamical 
Parameters 
Satellite  1 
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C.  Multisatellite  Covariance  Upper  Triangular  Halt  (Continued) 


Word  Number 

Format 

Definition 

• 

\ 

Dynamical 

• 

• 

• 

Parameters 

• 

• 

• 

Satellite  2 

24 

E 

Cov  XI,  T3  \ 

(N-  1)12+  1 

E 

Cov  XI.  Xn  j 

Dynamical 

(N-  1)12  +  2 

E 

Cov  XI,  Yn  ' 

,  Parameters 

• 

• 

Satellite  N 

# 

) 

(N)12 

E 

Cov  XI.  T3 

E 

Cov  XI,  TB,  1 

Time  Satelliet  1 

E 

Cov  XI.  TB, 

. 

E 

Cov  XI.  TB,  | 

• 

Cov  XI,  TB, 

• 

Cov  XI,  TB’  1 

Time  Satellite  2 

E 

Cov  XI,  TB~  ] 

• 

Cov  XI.  TBN  ' 

• 

• 

1  Time  Satellite  N 

E 

Cov  XI,  TBn 

| 

E 

Cov  X),  TB,  i 

I 

E 

Cov  XI,  TB, 

l  Time  Station  1 

• 

Cov  XI,  TB,  ; 

I 

• 

Cov  XI,  TBM 

• 

• 

Time  Station  M 

• 

•  ..  1 

(Master  Station) 

E 

Cov  XJ,TBM  i 

Row  2 

1 

E 

Cov  Y1  Y1 

2 

E 

Cov  Y1  Z1 

3 

E 

Cov  Y1  XI 

etc. 

etc. 

Row  3 

1 

E 

Cov  Z1  Z1 

2 

E 

Cov  Z1  XI 

3 

E 

Cov  Z1  Y1 

etc. 

etc. 
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to 
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APPENDIX  G 
ABCD  FILE 
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A  BCD  FILE 


This  file  contains  the  ABCD  rotation  matrix  from  CELEST  and  the  Greenwich  Hour 
Angles  (arguments  of  the  B  matrix).  The  first  record  is  a  header  and  the  remaining  records 
correspond  to  each  time  line  in  the  integration. 


RECORD  TYPE  1 


Word  Number  Type 


1  I 

2  F 

3  F 

4  1 

5  I 

6  I 

RECORD  TYPE  2 

Word  Number  Type 


1  I 

2  F 

3-11  F 

12  Double 


Precision 


Description 


Block  Number  (=  1 ) 

Epoch 
Step  Size 

KGC(3)  Not  Used  by  SAMSAP  Program 


Description 


Block  Number  (=  2,  3,  4  •••) 

Time  in  Seconds  from  Epoch 
ABCD  Transformation  Matrix  (Column 
Wise  3X3) 

Greenwich  Hour  Angle 
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APPENDIX  H 

DIAGNOSTIC  INFORMATION  FILE 


H-l 


DIAGNOSTIC  INFORMATION  FILE 


The  format  of  this  file  is  similar  to  that  of  the  Diagnostic  Information  File  (Tape  73) 
of  CELEST  (Appendix  C).  It  has,  however,  been  expanded  to  handle  the  multiple  satellite 
case  and  the  clock  parameters  required  by  the  SATRACK  program.  A  header  record 
containing  the  number  of  satellites  and  stations  has  been  added.  Record  types  2,  3,  4,  and 
5  appear  for  each  satellite  separated  by  an  end-of-file.  Record  type  4  appears  for  each 
output  time  for  each  satellite. 

RECORD  TYPE  1 


Word  Number 

Type 

Description 

1 

I 

Number  of  Satellites  (N) 

2 

I 

Number  of  Stations  (M) 

RECORD  TYPE  2 

Word  Number 

Type 

Description 

1 

F 

Year  of  Epoch 

2 

F 

Day  of  Epoch 

3 

F 

Seconds  of  Epoch 

4 

l 

Satellite  Number 

5 

I 

Scale  Factor  (=1) 

6 

I 

Rev  Number  (=  1 ) 

7 

F 

Year  of  Epoch 

8 

F 

Day  of  Epoch 

9 

F 

Seconds  of  Epoch 

10 

F 

Output  Delta  T  (step  size) 

11 

I 

Dummy  (=1) 

RECORD  TYPE  3 

A 

Word  Number 

Type 

Description  J / 

1 

F 

Last  Time  for  Which  Record 
Type  3  is  Written  (stop  time 
for  outputting  covariance 
matrices) 
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RECORD  TYPE  4 


Word  Number 

Type 

Description 

1 

F 

Time  in  Seconds  From  Epoch 

2 

F 

Sigma  Radial  Component  Satellite  1 

3 

F 

Sigma  Tangential  Component 

Satellite  1 

4 

F 

Sigma  Out-of-Plane  Component 

Satellite  1 

5 

F 

Sigma  Clock  Bias  Satellite  1 

6 

F 

Sigma  Clock  Drift  Satellite  1 

7 

F 

Sigma  Clock  Aging  Satellite  1 

8 

F 

Sigma  Clock  Bias  Station  1 

9 

F 

Sigma  Clock  Drift  Station  1 

10 

J 

■ 

F 

Sigma  Clock  Aging  Station  1 

J 

J 

3(M-  1)  +  7  +  1 

F 

Sigma  Clock  Bias  Station  M 

3(M-  l)  +  7  +  2 

F 

Sigma  Clock  Drift  Station  M 

3(M  -  l)  +  7  +  3 

F 

Sigma  Clock  Aging  Station  M 

Record  Type  4  is  output  for  each  time 
Record  Type  4  (for  each  satellite)  is  zero  filled. 

in  covariance  propagation  span.  The 

RECORD  TYPE  S 

Word  Number 

Type 

Description 

1 

2 

F 

F 

Not  Used  -  Should  Contain 

3 

F 

Dummy  Values 

An  end-of-file  is  written  after  Record  Type  5.  Record  Types  2,  3,  4,  and  5  are  then 
repeated  for  the  next  satellite.  A  double  end-of-file  is  written  after  Record  Type  5  for  the 
last  satellite. 


